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Abstract. Let D be a simply connected, smooth enough domain of R^. For L > consider the 
continuous time, zero-temperature heat bath dynamics for the nearest-neighbor Ising model on 

with initial condition such that a^: = — 1 if s G LV and = +1 otherwise. It is conjectured 
[24] that, in the diffusive limit where space is rescaled by L, time by and L — >■ oo, the 
boundary of the droplet of "— " spins follows a deterministic anisotropic curve-shortening flow, 
where the normal velocity at a point of its boundary is given by the local curvature times an 
explicit function of the local slope. The behavior should be similar at finite temperature T < Tc, 
with a different temperature-dependent anisotropy function. 

We prove this conjecture (at zero temperature) when D is convex. Existence and regularity 
of the solution of the deterministic curve-shortening flow is not obvious a priori and is part of 
our result. To our knowledge, this is the first proof of mean curvature- type droplet shrinking 
for a model with genuine microscopic dynamics. 
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1. Introduction 

Consider a thermodynamic system with two coexisting phases and imagine to prepare it 
in an initial condition where a droplet of one phase is immersed in the other phase. If the 
system undergoes a dynamics that does not conserve the order parameter, it is well understood 
phenomenologically [20] that the droplet will shrink in order to decrease its surface tension 
until it eventually disappears, and that (roughly speaking) the normal speed at a point of its 
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boundary will be proportional to the local mean curvature. Deriving such behavior from first 
principles, i.e. from a microscopic model undergoing a local (stochastic) dynamics, is a much 
harder task and this program was started only rather recently [24]. More precisely, what one 
expects is that, if the initial droplet is of diameter L, it will "disappear" within a time of order 
LP' (this behavior is sometimes referred to as "Lifshitz law"). Moreover, in the "diffusive limit" 
where L — )• oo and at the same time space is rescaled by L (so that the initial droplet is of 
size 0(1)) and time is accelerated by L^, the droplet evolution should become deterministic and 
follow some anisotropic version of a mean curvature flow. Anisotropy (i.e. the fact that the 
normal velocity will also depend on the local orientation of the droplet boundary) is expected 
when the underlying model is defined on a lattice, as will be the case for us. 

Up to now, mathematical progress on this issue has been rather modest, the main difficulty 
being that it is not clear how to implement the idea that the fast modes related to relaxation 
inside the two pure phases should decouple from slow modes related to the interface motion, 
which are responsible for the diffusive time scaling L^. 

A fairly well understood situation is that where the interface can be described by a height 
function and the bulk structure of the two phases is disregarded: this is possible (by definition) 
for the so-called "effective interface models" or Ginzburg-Landau interface models: for mod- 
els with continuous heights and strictly convex potential undergoing a Langevin-type dynamics, 
Funaki and Spohn [10] derived the full mean-curvature motion in the diffusive scaling. Another 
well-studied case is that of models with Kac-type potentials: in this case, mean-curvature mo- 
tion can be proven to emerge [5, 6, 15] in a limit where interaction range is taken to infinity at 
some stage, but in this limit there is no sharp interface separating the phases and the system 
becomes very close to mean-field. 

As for true lattice models, results are much more scarce. For instance, for the two-dimensional 
nearest-neighbor Ising model below the critical temperature, the best known upper bound on 
the "disappearance time" for a droplet of "— phase" immersed in the "+ phase" is of the order 
2^c(T) log L j22]j very far from the expected scaling. Recently, a weak version of the Lifshitz 
law was proven for the three-dimensional Ising model at zero temperature: the disappearance 
time of a "— " droplet is of order (upper and lower bounds), up to multiplicative logarithmic 
(in L) corrections [2]. When the dimension is higher than three (always at zero temperature), 
an upper hound for the disappearance time of order L^(logL)'^, for some constant c, was proven 
in [18]. 

In this work, we concentrate on the two-dimensional nearest-neighbor Ising model on the 
infinite square lattice. The dynamics takes a very simple form: each spin is updated with rate 
one and after the update it takes the same value of the majority of its neighbors, or the value 
±1 with equal probabilities if exactly two neighbors are "-I-" and two are "— ". In this case, 
the disappearance time of a large "— " droplet should be asymptotically given by one half its 
volume (number of "— " spins). Moreover, in the diffusive scaling limit the droplet boundary 
should be given by a deterministic curve 7(i) whose normal speed is given by the local (signed) 
curvature, times a function a{9) where 9 is the angle of the local normal vector. The function 
a(-) is explicitly known, see (2.4). In this two-dimensional setting, it is more natural to refer to 
such ffow as "(anisotropic) curve-shortening flow" (rather than "mean curvature flow"). 

Our main result (Theorem 2.2) is a proof of the curve-shortening conjecture (and, as a byprod- 
uct, of the Lifshitz law) when the initial droplet is convex. 

There are some previous partial results available on this problem. In [23] (cf. also [21]), Rost 
described (with the language of exclusion processes rather than spin systems) the scaling limit 
of the evolution when initially spins are "— " in the first quadrant of T? (infinite corner) and 
"-|-" elsewhere. This is a simple situation because the interface motion is mapped to symmetric 
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simple exclusion and is described by the associated height function at all times. In [24], Spohn 
described the scaling limit of the interface motion in a situation that more or less corresponds to 
the zero-temperature Ising model in an infinite vertical cylinder, with an initial condition such 
that the interface separating "+" from "— " spins can be globally described by a height function 
at all times (in particular, this cannot describe a droplet, and implicitly he has to modify the 
dynamics to guarantee that droplets do not appear later in the evolution). In [4], Chayes et al. 
proved the Lifshitz law (but not the curve-shortening conjecture) for a modified dynamics where 
updates which break the droplet into several droplets are forbidden. In [3], Cerf and Louhichi 
computed the "drift at time 0" of the droplet (for the non- modified dynamics), but their result 
does not allow to get information on the evolution for finite time t > 0. 

An important building block of our proof of the anisotropic curve-shortening conjecture is that, 
as was well understood by Spohn in [24], locally the interface can be (roughly speaking) described 
by the hydrodynamic limit of a certain zero-range process at the points where the tangent to 
the boundary is horizontal or vertical, and by the hydrodynamic limit of the symmetric simple 
exclusion process elsewhere. However, such correspondence is not exact due to updates that 
split the droplet into more than one connected components (see for instance Figure 4). In other 
words, the interface is not (even locally) the graph of a function. Also, it is a non-trivial task 
to patch together the various pieces of "local analysis" to control globally the evolution of the 
droplet. Both problems will be tamed by a sequence of monotonicity arguments, which are 
allowed because the dynamics conserves the stochastic ordering among configurations. 

In order to prove Theorem 2.2, we also need to know that a classical solution to the anisotropic 
curve-shortening flow exists up to the time where the droplet disappears, and (crucially) that 
such solution is sufficiently regular in space and time (i.e. that the curvature is a Lipschitz 
function of the angle and a continous function of time). To our surprise, we found that the 
existing literature on curve-shortening fiows does not provide global (in time) results for the fiow 
associated to the zero-temperature 2d Ising model. The reason is that the anisotropy function 
a(-) is not smooth (its derivative has jumps, reflecting the singularities of the surface tension 
at zero temperature), while the existing results assume that a(-) is at least C^, cf. [11, 12]. To 
prove existence, uniqueness and regularity of the solutions (cf. Theorem 2.1), we will regularize 
the function a(-) and then analyze the regularized fiow following the ideas of [11, 12]. Of course, 
it will be crucial to guarantee that all the estimates we need are uniform in the regularization 
parameter, which tends to zero in the end. 

The case where the initial "— " droplet is non-convex will be considered in future work. The 
additional difficulties are two-fold. First of all, from the analytic point of view, available global 
existence and regularity results for the solution of curve-shortening fiows with non-convex initial 
condition seem to be limited to the isotropic case where a(-) = 1 [14]. Secondly, due to the 
fact that the droplet will move at the same time outwards and inwards at different locations 
according to the sign of the curvature, various monotonicity arguments we use in the following 
will not work. 

2. Model and results 

2.1. Glauber dynamics and expected limiting evolution. Set Z* := Z + ^ := {x + 

(1/2) \x € Z}. We consider the zero-temperature stochastic Ising model on (Z*)^ with its 
usual lattice structure (x and y are linked if [a; — y] = 1 for the li norm). Is is a continuous 
time Markov chain {a{t))t>Q on the space of spin configurations on (Z*)^, 0, := {—1,1}^^ ^ . 
We write a{t) = (o"x-(t))a;e(z*)2 and for simplicity we write (Tx = — (resp. ax = +) instead of 
ax = -1 (resp. ax = +1). 



4 



H. LACOIN, F. SIMENHAUS, AND F. L. TONINELLI 



The transition rules are the following : for each site x G (Z*)^, the value Gx of the spin at x 
is updated independently with rate one. When the spin at a site is updated, it takes the same 
value as the spin of the majority of its neighbors, or the values ±1 with equal probabilities 1/2 if 
two neighbors have "+" spins and the other two "— " spins. That these rules yield a well-defined 
Markov chain even in infinite volume is a standard fact (cf. [21]). In the following (cf. (2.1)), 
we will consider only initial conditions where the number of "— " spins is finite. It is easy to 
realize that the spins outside the smallest square containing all the initial "— " spins stay "+" 
forever, so that in reality we have a dynamics on a finite volume and the question of existence 
of the process is trivial. 

We are interested in the evolution of the set of "— " spins for this Markov chain when the 
initial condition (t(0) is a large droplet, i.e. a finite connected set of "— " spins surrounded by 
"+" spins. In that case, almost surely, after a finite time r+, all the "— " spins have turned 
to "+" and the dynamics will stay forever in the all "+" configuration (which is an absorbing 
state). Our aim is to describe the evolution of the shape of the rescaled minus droplet on a 
proper (diffusive) time-scale. In the next paragraph we make that aim more precise. 

We consider a compact, simply connected subset V C [—1,1]^ whose boundary is a closed 
smooth curve. Given L S N we consider the Markov chain described above with initial condition 

^rm /-^ ifxG(z*)2nLP, 

c^x(O) = < ^ . (2.1) 

I +1 otherwise. 

In order to see a set of "— " spins as a subset of M^, each vertex x G (Z*)^ may be identified 
with the closed square of side- length one centered at x, 

:=x+[-l/2,l/2]2. (2.2) 

One defines 

AL{t):= U (2.3) 

{x: <T^(t) = -l} 

which is the "— droplet" at time t for the dynamics. The boundary of Ahit) is a union of edges 
of 1? (this is the only reason why we defined the Ising model on 1?). 

What was conjectured by Lifshitz [20] on heuristic grounds for the low temperature Ising 
model is that Ahit) should follow an anisotropic curve shortening motion: after rescaling space 
by L and time by and letting L tend to infinity, the motion of the interface between AL{t) and 
its complement should be deterministic and the local drift of the interface should be proportional 
to the curvature, with an anisotropic correction to reflect anisotropy of the underlying lattice. 
More precisely, one can formulate this conjecture as follows [24]: Let 7(t, L) denote the boundary 
of the (random) set (l/-L)^i(L^t). Then, for L — ?• oo, j{t,L) should converge to a deterministic 
curve 7(t) and the evolution of (7(i))t>o should be such that the normal velocity at a point 
X G 7(t) is given by the curvature at x, times an anisotropic factor a{9x), where 6x is the slope 
of the outwards directed normal to j{t) at x. The velocity is directed inwards at points where 
7(t) is convex and outwards at points where it is concave. The function a(-) should have the 
explicit expression 

""^^^ 2(|cos(^)| + |sin(0)|)2- (^-^^ 
In particular, the curve 7(t) should shrink to a point in a finite time 

Area{T>) _ Area{T>) 

J^- a{e)de ~ 2 



*0 - r2n 
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Note that the function a(-) is symmetric around and is periodic with period 7r/2, which 
reflects the discrete symmetries of the lattice Z^. It is important to note for the following that 
a(-) is C°° except at ^ = j7r/2, j = 0, ... ,3 where it is only continuous and its first derivative 
has a jump: indeed, a{9) ~ 1/2 — \0 — i7r/2| for 9 close to Z7r/2,z = 0, ... ,3. 

2.2. Results. 

2.2.1. Convex initial droplet. We prove the anisotropic curve shortening conjecture in the case 
where the initial droplet is convex and suitably smooth. Given a strictly convex smooth curve 
7 = &D in M^, we parameterize it following a standard convention of convex geometry (cf. e.g. 
[12] and Figure 1). For 9 G [0,27r] let v{9) be the unit vector forming an anticlockwise angle 9 
with the horizontal axis and let 

h{9) =su^{x-v{9),x e-i] (2.5) 
with • the usual scalar product in M^. 




Figure 1. A graphic description of the support function h. Given 0, consider 
the point x{9) of 7 that maximizes x ■ v{9) (it is unique if the curve is strictly 
convex). Then h{9) = x{9) ■ v{9), and k{9) is the norm of the curvature vector of 
7 (bold vector) at x{9). If the tangent to 7 at x exists it is normal to v{9) and 
\h{9)\ is the distance between the tangent and the origin. 



The function 9 1— >• h{9) (called "the support function") uniquely determines 7: 

V = no<e<2Ax G : X • v{9) < h{9)}. (2.6) 
With this parameterization, the anisotropic curve shortening evolution reads 

dthi9,t) = -ai9)ki9,t) 

h{9,0) = h{9) ^ ■ ' 

where, for a convex curve 7, k{9) > is the curvature at the point x{9) G 7 where the outward 
normal an anticlockwise angle 9 with the horizontal axis and the t-derivative is taken at constant 
9 (see [12, Lemma 2.1] for a proof of (2.7)). Of course h{-) is the support function of dV. 

In general, even proving the existence of a solution of (2.7) with a(-) given in (2.4) is non- 
trivial, since a(-) has points of non-differentiability and the existing literature (e.g. [11, 12]) 
usually assumes that a(-) is at least C^. 
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Our first result is 

Theorem 2.1. Let V C [-1,1]^ he strictly convex and assume that its boundary 7 = &D is 
a curve whose curvature [0, 27r] 3 6 k{9) defines a positive, 2'K-periodic, Lipschitz function. 
Then there exists a unique flow of convex curves {'y{t))t with curvature defined everywhere, such 
that 7(0) = 7 and that the corresponding support function h{6,t) solves (2.7) for t > and 
satisfies the correct initial condition h{9,0) = h{9). The curve 7(t) shrinks to a point xj G 
at time tf = Area{T>)/2. For t < tf, 7(t) is a smooth curve in the following sense: its curvature 
function k{-,t) is Lipschitz and bounded away from and infinity on any compact subset of 
[0,tf). 

We let P(t) denote the convex closed set enclosed by j{t) (of course, V{t = 0) = V). Also, 
we use the convention that T>{t) = {xj} if t > tj. 

For 5 > let B{x,6) denote the ball of radius 6 centered at x and for any compact set C cM? 
define 

:= U Bix,6), := B{x,6)\ . (2.8) 

Note that P(t)(^) = 5(x/, b) and P(t)(-^) = if t > t/. 

An event Bi^ is said to occur with high probability (w.h.p.) if lim/,_j.oo P{Bl) = 1. 

Theorem 2.2. Under the same assumptions on T> as in Theorem 2.1, for any 5 > one has 
w.h.p. 

V^-^\t) C ^AiiL^t) C V^^\t) for every 0<t<tf + 6 (2.9) 

AiiL^t) = for every t>tf + 5. (2.10) 
In particular, one has the following convergence in probability: 

lim ^ , , = -. (2.11) 

L^oc L^Area{V) 2 ^ ^ 

The reason why in Theorems 2.1 and 2.2 we do not content ourselves with, say, initial C°° 
curves is that, as we see in next section, there is a very natural initial condition whose curvature 
function is only Lipschitz and not (and stays so at later times). 

Theorem 2.2 does not apply directly if one considers V = [0, 1]^ or any other non-smooth 
convex set. However, approximating T> from above and below by smooth compact sets and 
using monotonicity (cf. Section 2.3), one sees easily that (2.11) holds in any case. In particular, 
the disappearance time of an L x L square droplet is with high probability L^/2(l + o(l)). 

Theorems 2.2 and 2.3 tell us that for our choices of initial configuration, the disappearance 
time of the minus droplet is non-random at first order. This implies that the variation distance 
of our Markov Chain from equilibrium (which is concentrated on the all-plus configuration) 
drops abruptly from 1 to around time L'^tj within a time- window of width o{L'^) <^ L^tf 
(we conjecture that the correct order of the window should be 0(L^/^)). This is a particular 
instance of a phenomenon called cut-off (cf. [7] and [19]). 

2.2.2. Scale-invariant droplet. A particular case of Theorem 2.2 is that where the initial con- 
dition is scale invariant, i.e. when the limiting evolution {'y{t))t is a homothetic contraction. 
Consider the function 



/o: 



1 1 



X 



3 X ^ fo{x) = ^ Uax e^"* dt - e^"^ ^ , (2.12) 







CURVE-SHORTENING EVOLUTION FOR THE 2D ISING MODEL 



7 



where a is the unique positive solution of 

4V2ae 

and 







(2.13) 



P = -\/2e-" < 0. (2.14) 

Note that /o is C°° , positive, concave, symmetric around and increasing on [~"^)0]- We 
denote by (ei, 62) the canonical basis of and (fi, f2) = ( ^^^^ ' ^'-J^^ ) the image of (ei, 62) by 
the rotation of angle — 7r/4. We also define the curve 71 to be the graph of /o in the coordinate 
system (fi, £2), i.e. 

71 := |xfi + /o(x)f2 I -i=,i= |. (2.15) 

If Si (resp. 5*2) denotes the symmetry with respect to the axis ei (resp. e2) one defines the 
closed curve 7 by 

7 = 71 U (Si7i) U (S271) U {{Si o ^2)71). (2.16) 
In the sequel ^ denotes the compact, convex set enclosed in 7, see Figure 2. 





1 














7 


/ V2 \ 


62 k{x)/ \ 










1 

|V2 



















Figure 2. The curve 7 = 33! and the coordinate systems (ei,e2) and (fi,f2). 



One can check that the curvature function 1— >• fe(0) of is Lipschitz and bounded away 
from zero, but not differentiable at = 111/2,1 = 0,1,2,3. In this case. Theorem 2.2 can be 
formulated as follows. 



Theorem 2.3. Assume that V 



K For any r] > 0, w.h.p, 



(Vl - 2at - r])& C -AL{tL^) C (Vl - 2at + r])Q for every t>0 (2.17) 
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where we work with the convention that y/x = for x < and that x& = for x < 0. Moreover, 
one has the following convergence in probability: 

lim ^, = a lim = -. (2.18) 

L^oo Area{L&) L^oo 2 

It is easy to check, using Lemma 3.6 below and a couple of integrations by parts, that 
Area{S!) = 1/a, yielding the first equality in (2.18). The expression (2.12) for the invariant 
shape appears also, although with different notations, in the recent work [17]. 

2.3. Graphical construction of the dynamics and monotonicity. Before starting the 
proofs, we wish to give a construction of the Markov process (called sometimes the graphical 
construction) that yields nice monotonicity properties. We consider a family of independent 
Poisson clock processes {t^)x(^7?- More precisely, to each site x G (Z*)^ one associates a random 
sequence (independently from other sites) of times (T^)n>0; that are such that Tq = and 
{'T'n+i ~ Tn)n>{) are IID exponential variables with mean one. One also defines random variables 
{Un,x)n>o^x G (Z*)^ that are IID Bernoulli variables of parameter 1/2, with values ±1. 

Then given an initial configuration ^ G { — 1, l}^^*-*^ one constructs the dynamics (J^{t) starting 
from (T^(0) = as follows 

• (<7x(i))t>o is constant on the intervals of the type [t^,t^+i). 

• (Tx{Tn) is chosen to be equal to ±1 if a strict majority of the neighbors of x satisfies 
cry{T^) = ±1, and Un,x otherwise (this definition makes sense as, almost surely, two 
neighbors will not update at the same time unless with zero probability.) 

This construction gives a simple way to define simultaneously the dynamics for all initial 
conditions (we denote by P the associated probability). Moreover this construction preserves 
the natural order on { — 1, +1}*^^*) , given by 

C>C'^^x>C'x for every x G (2.19) 

(this order is just the opposite of the inclusion order for the set of "— " spins, which is therefore 
also preserved). Indeed, if ^ > C with the above construction, one has P-a.s. 

yt>0 a^it)>a^'{t). (2.20) 



3. Local interface dynamics 

One problem one has to deal with when proving mean curvature motion for the whole droplet 
is that even though initially the interface between "+" and "— " (i.e. the geometric boundary of 
the set Al{0)) is a simple curve, it can later split to form several loops. In fact, as a byproduct 
of our results, we will obtain that, with large probability, only very small extra loops can be 
created. We will tackle this problem by introducing some auxiliary dynamics that do not allow 
creation of new loops and stochastically compare to the original one. 

A second problem is that the interface that one has to control is not exactly the graph of 
function, for which it would be easier to describe the macroscopic motion using partial differential 
equations. We begin by studying two dynamics for which the interface is indeed a graph, and 
which have locally the same large-scale behavior as the true evolution. It is more natural to 
introduce these dynamics as dynamics on interfaces rather than dynamics on spins. Our work 
then will consist of glueing together the "local results" of Theorems 3.2 and 3.4 to get Theorems 
2.2 and 2.3. 
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3.1. Local dynamics away from the poles and simple exclusion process. The first 
auxiliary dynamics is used to control the evolution of the boundary of {l/L)ALitL'^) away from 
the points (the poles) where the tangent to the deterministic curve 7(t) is either horizontal or 
vertical. The evolution near the poles will be analyzed via a second auxiliary dynamics, see 
Section 3.2. 

Given two positive natural numbers M, N consider the state-space 0,m,n of nearest-neighbor 
directed paths of length L := M + N with M steps up and N steps down: 

^M,N = {(/ix)x-e{o,...,Af+7V} G I \hx+i -hx\ = 1, /lo = 0; Hm+n = M - N}. (3.1) 

Given h G i^M,N and x G {1, . . . ,L — 1}, we denote by h^^^ the path with a corner "flipped" 
at X defined by h^^'^ = hy for all y ^ x and 

hx — 2 ii hx±i = hx — 1, 
4"^ ■.= lhx + 2if hx±i = hx + l, (3.2) 



hx if \hx+i - h. 



X I 



The dynamics ^m,n we consider is the one that flips every corner with rate 1/2. More 
precisely it is the Markov chain whose generator C is defined as 



1 ^"^ 

Cf{h):=-Y,U{h'^''^)-f{h)). yf:nM,N^^. (3.3) 

x=l 

We denote by {h{t))t>o the trajectory of the Markov chain started from initial condition 
h{0) := h^ G ^M^N. 

Remark 3.1. Note that this dynamics is in one-to-one correspondence with the Ising dynamics 
on a rectangle N x M with boundary conditions on two adjacent sides and "— " boundary 
conditions on the two opposite sides, provided that the initial configuration is such that the 
length of the — /+ boundary is Af -|- (i.e. the minimal possible length). More precisely (see 
Figure 3) the correspondence is obtained by taking the graph of h, rotating it by 7r/4 and 
rescaling space by a factor (so that squares have side- length one on the right picture). 
Note that we are implicitly identifying an element h G ^Im,n with a continuous function F : 
[0,M + N] i-> M such that F{x) = h^ fov x = 0,1, . . . , M + N and F(-) is linear on intervals 
{n,n + 1) with integer n. 

This corner-flip dynamics has been widely studied (see e.g. [25]) and can be mapped to the 
symmetric simple exclusion process (SSEP) on a finite interval (just say that there is a particle 
at X = 0, . . . , M + N — 1 if and only if hx+i — hx = +1, and check that dynamics in terms of 
particles coincides with that of SSEP). From hydrodynamic-limit results, it is quite clear that 
the rescaled version of h when M, N tends to infinity should satisfy the heat equation (see [16, 
Section 4.2.] for an account on hydrodynamic equations for the exclusion process). However, 
we have not found in the literature a proof of the following precise statement we need (we give 
a concise proof of it in Section 7): 

Theorem 3.2. Given a 1-Lipschitz function (fP : [0,1] i— )• M with 4P{0) = 0, let {h{t))t>o the 
dynamics starting from initial condition h^ G Qml,Nl given by 

hi := 2[L/(x/L)/2J for even x 

hi := 2[(L/(a;/L) - l)/2j + 1 for odd x ^^'"^^ 
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M - N 




M + N 



Figure 3. One-to-one correspondence between the dynamics in a rectangle with 
mixed boundary conditions and the corner-flip dynamics on paths. A possible 
spin update together with the equivalent corner-flip are represented 



where (j) : [0, 1] x 



(M^ and Nl are implicitly fixed by L and ). For all T > and e > 0, w.h.p. 

sup y - L(l){x,t)\ < e (3.5) 

te[0,T],a;e[0,l] ^ 

[ is the solution of the Cauchy problem 
dt^{x, t) = \dl(t){x, t) \ft > 0, Vx G (0, 1), 

(t){0,t) = 0, (^(l,t) = /(l) Vt>0, (3.6) 
^(j){x,0) = (j)^{x) VxG(0, 1). 

3.2. Local dynamics around the poles and a zero-range process. For the definition of 
the second auxiliary dynamics, we use the same notation as in the previous section, but no 
confusion should arise as the proofs will be given in two independent sections. The state space 
is 

nL:={h:{-L,...,L+l}^Z}. (3.7) 
For h G 0,L and x G {—L + 1, . . . ,L} define /i"*"'^ (resp. h~'^) as the configuration such that 
hy'^ = hy ii y ^ X and hx'^ = + 1 (resp. hx'^ = hx — 1). We consider the Markov chain 
{h{t))t>o started from some ^ Ql and with generator C defined by 

1 ^ 

^/w = o E c+'^{hmh^'n-fm+c-'^h){f{h-^^)-f{h)) m 



-L+l 



where 

c ' 

c" ' [n) = ±{/i^+i</i4 -i-{h,_i</i4- 

Note that the values H-l and are fixed in time and should be considered as boundary 

conditions. 



(3.9) 



Remark 3.3. This dynamics corresponds to the motion of the interface for a modified Ising 
dynamics in a vertical strip of width 2L with the following boundary condition: spins on the 
left (resp. right) boundary of the system are "+" if and only if their vertical coordinate is 
larger than (resp. /il+i)- The dynamics is modified in the sense that updates are discarded 
if after the update the boundary between the "— " and "+" domain is not a simple (open) 
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curve (see Figure 4). It is at times more convenient to identify h G with a cadlag function 
H : [—L — l/2,L + 3/2] i— >• Z which equals identicahy hn on intervals [n — 1/2, n + 1/2) for integer 
n. 



+ spins 



spins 




Figure 4. An example of spin update that splits the interface into two discon- 
nected components. The interface dynamics presented in this section does not 
allow this kind of move. 



Another way to interpret this dynamics [24] is to look at the gradients rjx = h^+i — h^: one 
recognizes then a zero-range process with two type of particles {\i r]x = n > Q we say there are 
n particles of type A at x, if rjx = —n < we say there are n type-B particles). Each particle 
performs a symmetric simple random walk with jump rate 1 /(2n) (with n the occupation number 
of the site where the particle sits) to either left or right and particles of different type annihilate 
instantaneously when they are at the same site. See Figure 5. 




Figure 5. Correspondence between interface dynamics and zero-range process. 
Arrows represent possible motions for the interface and their representation in 
terms of particle moves. When an A particle jumps on a B particle (green arrow) 
both annihilate. 



In [24, Appendix A], this dynamics was considered but in a periodized setup. A scaling 
limit result was given but the proof there is somewhat sketchy. Here we adapt the proof to the 
non-periodic case and write it in full details. 

Consider (fp : [—1, 1] i— t- R a function with ^''(1) = = 0. We further assume that 

(j)^ has a finite number of changes of monotonicity. Define <l>o : {—L, . . . , L + 1} i— > M as 

^o{x) := L^^{x/L) (3.10) 

and /i° : {-L, . . . ,L + I] ^ Zhy 

hl:=[^o{x)\. (3.11) 
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We define $ : {— . . . , L + 1} x M+ — M as the solution of the following Cauchy problem: 

'dt^{x,t) = \ \a{q^{t)) - a{q^.xiM , 

< $(L + l,t) =$(-L,i)=0, (3.12) 
$(x,0) =^o(a;) 
for every t > and x G {— -i^, . . . , L + 1}, where (y{u) = u/(l + |n|) and 

q^{t): = ^{x^\,t)-^{x,t). (3.13) 

The result we state now is slightly weaker than Theorem 3.2 as it allows to control the profile 
h only at a fixed time and not on a whole time interval. 

Theorem 3.4. Given (fP as above, consider (/j(t))t>o the dynamics described by (3.8) with initial 
condition as in (3.11). Then for any t, the following convergence holds in probability 

lim max \\h^{LH) - ^{x,L^t)\ = (3.14) 
L^oox&{-L,...,L+l} L ^ 

It is quite intuitive that one should have that -j;$([-LxJ , L^t) — )• (j){x,t) for any x G [—1,1], 
where (f) : [—1,1] x — )• M is the solution of 

ut(pyx,i) - 2 

0(1, t) =0(-l,t) = O (3.15) 
[0(x,O) =/(x) 

for t > and x G (—1,1). The particular form of the non-linearity of this PDE makes the 
convergence question non-trivial, but fortunately Theorem 3.4 together with a comparison with 
the heat equation (cf. Section A. 5) turns out to be sufficient for our purposes. Indeed, define 
4) : [—1, 1] X ]R+ — )• M to be the solution of 

^dt^{x,t) =\dl4){x,t) 
4>{l,t) =0(-l,t) = O (3.16) 
>(x,0) =ctP{x). 

Then 

Corollary 3.5. Let (fp be as above, and assume further that it is concave with JJc^x'^'^lloo < 
For every t > and every e > the following inequality holds w.h.p. 

4>ix/L, t)-e< jh^iL^t) < 4>{x/L, (1 + r/)"^*) + e for every x G {-L, ...,L + 1}. (3.17) 

Proof. The result follows by combining Theorem 3.4, Proposition A. 9, and by taking limits of 
rescaled versions of <I>i and <I>2 in (A. 55) when L tends to infinity (cf. Lemma 7.1). □ 

3.3. About the scale-invariant shape. Now that we know how the interface should evolve 
locally (from Theorems 3.2 and 3.4) it is possible to explain why ^ should be scale invariant. By 
symmetries of the problem and the fact that motion is driven by curvature, the scale-invariant 
shape should be convex symmetric around the axes Mei, Me2. Therefore it is enough to consider 
the boundary of the intersection of & with the first quadrant. 

From Theorem 3.2, if / is a Lipschitz function and d& is the graph of / in the coordinate 
system (fi,f2), the initial drift in the f2 direction is (1/4)5^/, where the factor 1/4 (instead 
of 1/2) is due to the fact that in the correspondence between Ising dynamics and dynamics of 
nearest-neighboring paths, space has to be rescaled by ^/2, cf. Remark 3.1). One the other 
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hand, the homothetic contraction of a shape D of initial velocity a gives an initial drift of the 
interface in the f2 direction 

a{-f + xdj). (3.18) 
That leads to the partial differential equation 

dlf = Aa{-f + xd,f). (3.19) 

Next we impose the correct boundary conditions on /: 

• We fix the scaling by imposing that the point (1,0) (and therefore also (0, 1), (—1,0), (0, — 1)) 
belongs to d^^. This gives 

/ (±1/^) = I/V2. (3.20) 

• To guarantee that the curvature of dS^ is well defined at the point (0, 1) we have to 
impose 

d^f (-I/V2) = -dj (l/V2^ = 1. (3.21) 

We finally notice that 

Lemma 3.6. The function Jq defined in (2.12) is the unique solution of the Cauchy problem 
(3.19)-(3.20)-(3.21) for X S (— 1/\/2, +l/-v/2). For other values of a the above problem has no 
solution. 

Proof. Uniqueness of the solution is standard from theory of ordinary differential equation. The 
rest is just a matter of checking. □ 

3.4. Organization of the paper. Instead of proving directly Theorem 2.2 and then deducing 
Theorem 2.3 as a corollary, we decided for pedagogical reasons to give first the proof in the 
case of the scale-invariant droplet and then to point out what needs to be modified in the more 
general case of a convex droplet. The reason is that, this way, we can easily separate the question 
of comparing the stochastic evolution with the deterministic one (which works more or less the 
same in the two cases but is simpler for the invariant droplet, due to its symmetries) from the 
analytic, PDE-type issues which appear only in the general case. 
The paper is therefore organized as follows: 

• in Section 4, we show that to prove Theorem 2.3 it is sufficient to have a good control 
on the continuity of the interface motion (Proposition 4.2) and a result on the evolution 
after an "infinitesimal time" eL^ (Proposition 4.1). Such crucial results are proven in 
Sections 4.3 and 4.4; 

• in Section 5 we first prove Theorem 2.1 on the existence of a solution to (2.7), and then 
we prove Theorem 2.2 via a suitable generalization of Propositions 4.2 and 4.1; 

• finally, the hydrodynamic limit results of Theorems 3.2 and 3.4 are proven in detail in 
Section 7 and the Appendix A respectively. 

4. Proof of Theorem 2.3: evolution of the scale-invariant droplet 

4.1. Reducing to an "infinitesimal" time interval. We decompose the proof of Theorem 
2.3 into two propositions. The first (and the main one) says that after a time eL^ the droplet 
looks very much the same but contracted by a factor (1 — ae + o{e)). 
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Proposition 4.1. For all 5 > there exists eo{6) > such that for all < e < eo{5), w.h.p., 

AiiL^e) C {l-e{a-6))L&, (4.1) 

and 

AiiL'^e) D {l-e{a + 6))L&. (4.2) 
The second proposition controls continuity in time of the rescaled motion: 
Proposition 4.2. For every 6 > 0, w.h.p., 

AiiL'^t) C (1 + 6)LS! for every t > 0. (4.3) 
Moreover, for every 5 > there exists e > such that w.h.p 

AiiL'^t) D (1 - 5)L& for every t £ [0, e]. (4.4) 

Proof of Theorem 2.3 assuming Propositions 4-1 and 4. 2. Given r] fix b small enough and e < 
eo(^)- Then using (4.1) one gets that w.h.p. 

AL{L^e) C (1 - (a - S)e)LQ. (4.5) 

Let (^^^^(-L^t))t>o denote the evolution of the set of "— " spins for the dynamics started from 
initial condition "— " on (1 — e(a — 5))LSi and "+" elsewhere. Then using the Markov prop- 
erty and monotonicity of the dynamics, one can couple the dynamics {Al{LP'{£ + t)))t>o and 
{A^l\L'^t))t>o such that on the event (4.5) 

AL{L^{e + t)) C A^l\L^t), for every t > 0. (4.6) 

Therefore, after conditioning to the event in (4.5) and using (4.1) for (1 — (a — 6)e)L instead of 
L, one gets that w.h.p.: 

AL{Lh{l + (1 - (a - 6)ef)) C A'^1\l'^{1 - (a - 6)efe) C (1 - (a - 5)efL&. (4.7) 

Using this argument repeatedly one gets that, w.h.p., for all k G [1, (here and in the sequel 

we assume e~^/^ to be in N) 

AULhk) C (1 - (a - 5)e)''L& (4.8) 

where tk is defined by 

fc-l , /, / .-N ^^2k 



. . ^V-/-| / .X x2i ^ l-(l-(a-(5)g ^ 

^ 1 - (1 - (a - d)eV 

'1=0 

One remarks that for all values oi k G [l,e~^/^]. 



(l-(a-.).)- = ,/l-^-(^-(^-("-^^-)^) 



(4.10) 



= ^l-2{a-6)tk + tkOie). 
As (tfc)fe ^ is bounded above, there exists C > such that for every fee [1,6-3/2] 



AiiL^tk) C I - 2(a - 5)tk + Ce)LS)C (VI - 2atk + ??/2)L^ (4.11) 



w.h.p. where the second inclusion holds provided that e and 5 are small enough. Combining 
(4.11), Proposition 4.2 and stochastic coupling, one gets w.h.p. that, for every k G [0, £-3/2] and 

Al{LH) C (VI - 2atfc + (37?/4))L^ C (Vl - 2at + r/)L^ (4.12) 
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and that, w.h.p., for every t > 

AiiLh) C (v/l-2ai,-3/2 + (37?/4))L^ C r]L^. (4.13) 

This ends the proof of the upper inclusion in (2.17) (note that approaches 1/(2q!) for e, 5 

small). Moreover, (4.13) and stochastic domination implies that, for some constant C, w.h.p. 

r+<^(l + Cr?2). (4.14) 

Indeed, it is known from [9] that a droplet of minus spins of linear size r/L disappears within a 
time r_|_ which w.h.p. is upper bounded by Cyf'Lp'. 

The lower inclusion in (2.17) and the lower bound on r+ are proved in an analogous way using 
(4.2) instead of (4.1) and (4.4) instead of (4.3). Note that using (4.4) we have to take care to 
choose e small enough but it is possible as tk — ifc-i is a non-increasing function of k and e. 

□ 

4.2. Strategy of the proof of Proposition 4.1. Our aim is to use Theorem 3.2 to control the 
motion of the interface away from the "poles" and Theorem 3.4 (or more precisely Corollary 3.5) 
to control the motion of the interface close to the "poles" . It is therefore crucial to compare the 
local SSEP or the zero-range dynamics introduced in Sections 3.1 and 3.2 to the true evolution 
of the boundary between "+" and "— " spins. 

As we have already discussed at the beginning of Section 3, however, there exists no exact 
mapping between the evolution of the height function associated to the two particle processes and 
the evolution of the +/— boundary, since the original "— " droplet can break into more droplets 
and, strictly speaking, the interface cannot be described, even locally, as a height function. The 
way out is that, thanks to monotonicity arguments and to the a priori "continuity" information 
provided by Proposition 4.2, we can remove certain updates of the Markov Chain, e.g. freeze 
certain spins to their initial value. This way, we can show that locally the interface can be 
stochastically compared to the height function associated to the SSEP (or to the zero-range 
process close to the poles). Of course, the detail of the "update removal procedure" is quite 
different according to whether we want to prove an upper or a lower bound on the "— domain" . 
For instance, if we want an upper bound we are allowed to freeze "— " spins or to change some 
"+" into "— " spins in the initial condition (this is fine thanks to monotonicity) and at the 
same time we can freeze the spins outside (1 + 5)LS' to "+" (this is not allowed directly by 
monotonicity, but (4.3) guarantees that such spins stay "+" for all times anyway, w.h.p.). If the 
"update removal procedure" is performed suitably, the effect is that the various portions of the 
+/— interface (away from and close to the poles) then become independent and evolve exactly 
like the height functions of the SSEP /zero-range process. 

The approach outlined here will be also used in Section 6 in the case with general convex 
initial condition (the generalization of Proposition 4.1 is Proposition 6.2). 

4.3. Upper Bound: proof of (4.1) and (4.3). The inclusion (4.1) can be rewritten in the 
following manner, which is more convenient for the proof: for any positive 6, for all £ small 
enough, w.h.p. 

a^{eL^) = + for every x G [(1 - e(a - 5))L^]' . (4.15) 

Given 5, we fix a value of ^ which is small enough (depending on 5 in a way that is specified 
in Section 4.3.2) and set (cf. Figure 6) 

M{e,0 ■.= {{x,y)eR', x ^ ^ ^^nd y ^ - e{a - 6)m 

o (4-lo) 
iV(e,C) := {(x,2/) G y ^ and - e ^ X ^ \ [(1 - e(a - 5))i^] . 
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Remark that for any e > 0, M,N and their successive images by rotation of angle 7r/2, vr, 3tt/2 
form an 8-piece cover of the complementary set [(1 — e{a — 5))!^Y. 




Figure 6. The light-colored (resp. dark-colored) zones correspond M(e,^) 
(resp. N(e,(^)) and its rotations. Together, they form a partition of the comple- 
ment of (1 — e(a — 6))& (white central region). 

As the dynamics and the initial shape are invariant under these same rotations, (4.15) is 
proved if we can show that for e small enough, w.h.p. 

a^{eL'^) = + for every x £ LM{e, £,) (4.17) 

o-^(eL2) = + for every x G LN{e,(). (4.18) 

The above new formulation of (4.1) is very convenient as it allows to consider separately the 
dynamics close to the poles and away from them. 

4.3.1. Proof of (4.17) . For any L > 0, we consider the dynamics which has initial condition with 
"— " spins in and "-|-" otherwise, and the same generator as the original dynamics except that 
spins on the sites inVi := {±^}x {—L + ^, ■ ■ ■ , L — i} and on V2 := {—L + ^, ••• ,L — i}x{ibi} 
are "frozen to — ". (The construction of the dynamics is the same as in Section 2.3, except that 

there is no update for these sites). We denote (c^^^(t))t>o the evolution of this dynamics and 
define 

A^lHt):= U C (4.19) 

{x, ai'\t) = ~l} 

The graphical construction of Section 2.3 gives a natural coupling of a and a^^^: 

Aiit) C A^l\t), for every t ^ 0. (4.20) 

The advantage of the "freezing procedure" is that then the evolution in the four quadrants 
of (Z*) becomes independent. The reason is that the spins on sites ({±i} X Z*) \ Vi and 
(Z* X {±^}) \ V2 a-re "+" for all times, since they always have at least three "-I-" neighbors. 
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The set A^l\t) nM.1 is a Young diagram (i.e. a collection of vertical columns of width 1 and 
non-negative integer heights, with heights non-increasing from left to right) for all t > and 

we can thus consider dA^^ (t) Ci as the graph of a (random) piecewise linear function in the 
coordinate system (fi,f2), that we denote by Fi{-,t). Equation (4.17) is thus proved (for any 
choice of ^) if one proves that for any u < 2"^/^, for any e small enough, w.h.p. 

FLix,eL^) < Lf{{a - 5)e,x/L) for every x £ {-pL.vL), (4.21) 

where f{-,t) is the function whose graph in the coordinate system (fi,f2) is given by the in- 
tersection of the boundary of (1 — t)^ with the half-plane {{x,y) G E?, {x,y) • f2 ^ 0} (the 
domain of definition of f{t, •) depends on t but includes [— 2~^/^, 2~^/'^] for t small enough). By 
definition of one has f{-,0) = /o(-) (recall the definition of /o in (2.12)). 

In practice, to prove (4.17) one has to prove (4.21) with u such that l/\/2 — = ^/\/2 + o(0 
for ^ small. The reason is that the point of dS^ with horizontal coordinate ^ and positive vertical 
coordinate (in the coordinate system (ei,e2)) has horizontal coordinate —(1 — ^)/\/2 + o(C) iii 
the (fi,f2) coordinate system. 

As explained in Remark 3.1 and on Figure 3, the function FL{-,t), up to space rescaling (by 
a factor ^/2) undergoes the corner fiip-dynamics of Theorem 3.2. Thus the scaling limit of 
satisfies the heat-equation or more precisely we have the following convergence in probability 
for every fixed T > 0: 

^FL{xL,tL'^) - g{x,t) 



lim sup sup 



L 



(4.22) 



where g is the solution for i > and x G (— 1/\/2, l/\/2) of 

'dtgix,t) = ldlg{x,t) 

5(-,0) = /o(-) (4.23) 
.<7(-^,i)=5(^,t) = ^. 

Note that the above result plus equation (4.20), plus the fact that g is decreasing in t (since 
it stays concave through time) gives (4.3) of Proposition 4.2 for every t < T < oo. Moreover, 
according to [9, Theorem 1.3], the disappearance time r+ of the minus-droplet is 0{L?') with 
high probability, so that (4.3) also holds for t > T provided that T was chosen large enough. As 
a byproduct, we have proven (4.3). 

Concerning (4.17), in order to prove (4.21) we are reduced to show that for every u G (0, l/\/2) 
and every x G (— z^, u), 

g{x,e) < fix,{a-6)e). (4.24) 

This is a consequence of the way /o was determined (see (3.19) and discussion in Section 3.3). 
First we notice that the time derivative of g is uniformly continuous away from the boundary 
points ibl/-v/2: 



Lemma 4.3. For any < v < 



limsup{\dtg{x,s) — dtg{x,0)\,s £ [0,t] and x £ = 0. (4.25) 

Proof of Lemma 4-3. This is well known but we sketch a probabilistic proof for the sake of 
completeness. Let (-Bt)i > o denote a standard Brownian motion starting at x G [— 2~^/^, 2~^/^] 
(with the associated expectation denoted by E^) and let T denote the hitting time of {ibl/-v/2}- 
One has 

dlg{x,t) = E, [a^/o(i?i)l|i<^|] . (4.26) 
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We can thus rewrite (4.25) as 

limsnp{\E^{dlfo{Bs) - d''joix))\,s G [0,t] and x G [-1^,1^]} = (4.27) 

and we can conclude using the uniform continuity of i9^/o on [— and well-known continuity 

properties of the Brownian motion. Remark that (4.25) would not hold with u = l/\/2 because 
of boundary effects: for t > one has that d^g{x,t) approaches zero as x approaches ibl/\/2, 
since Px{T > t) ^ when X ^ ±l/^/2. □ 

Therefore, for every r/ > arbitrarily small we have for all x in (— z^, i^), if e is small enough, 

g{x, e) < fo{x) + E {^^dlUx) + r?) . (4.28) 

We are left with proving that for e small enough and x in {—v, v) 

Ux) + e {^^dlUx) + r?) < /(x, (a - 5)e). (4.29) 

From the definition of /(•, t) as graph in (fi, £2) of the boundary of (1 — t)^, we get that if e is 
small enough, uniformly for all x G {—u,v), 



fix, (a - 5)e) = [!-(«- 6)e]f ^-^—^—^,oj (4.30) 

= foix) + {a- 6)e {xdJo{x) - fo{x)) + 0{e^). (4.31) 
Now recall that /o satisfies equation (3.19), so we are reduced to check that for all x G {—i^, v), 

V + S{xd,fo{x) - foix)) =V + ^dlfo < (4.32) 

which holds provided r/ = r]{6) is small, since d'^fo{-) is negative and uniformly bounded away 
from zero. Equation (4.17) is proven. □ 

4.3.2. Proof of (4.18). The method is similar to the one we used for (4.17), the main difference 
being that, via a chain of monotonicity arguments, we analyze the evolution of the portion 
of interface near the "poles" by comparing it to the interface dynamics of Section 3.2 (which 
coincides with the height function of the zero-range process with two types of particles) instead 
of the "corner-flip dynamics" . 

Denote by h{-,t) the function whose graph in the coordinates system (ei,e2) is given by the 
intersection of (1 — t)d!^ with the upper half-plane ffi x M_|_, and /io(") = h{-,0). Note that /iq 
is C°° on ( — 1,0) and on (0,+l) by the definition of ^. The boundary condition (3.21) ensures 
continuity of the first derivative of /iq at zero {dxho{0) = 0); the reader can check that ho has 
also continuous second derivative and that 

dlhoiO) = -i=9^/o(-l/V^) = -2a, (4.33) 

but that the third derivative exhibits a discontinuity in 0. 

Recall that ^ is the positive constant appearing in (4.16). We set h : [—4^,4^] 1— )• M to be the 
function defined by the following conditions: h = ho on [—2^, 2^], h is affine on [—4^, —2^] and 
on [2^,4^] and the derivative dxh{-) is continuous on (—4^,4^). Since /lo(') is strictly convex, 
we have ho{x) < h{x) with strict inequality outside [—2^,2^]. Define also the following subsets 
of M2 (cf. Figure 7): 

:= [4^, 00) X ^{Ai), 00) J2 := (-00, -4^] x [/i(4e), 00). (4.34) 
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To avoid notational complications with integer parts, we assume that L/i(4^) and 4L^ belong to 




Figure 7. In the set L{J^yj ,P) spins are frozen to "+" while in the dashed 
region they are is frozen to "— " . The initial condition is "+" in the dark-colored 
region and "— " in the light-colored one. The boundary separating dark/light 
region is determined by the function h{-). 



First of all, observe that, thanks to (4.1), we can freeze the spins in L{J^ U J^) to their initial 
value "-|-" and, w.h.p., the dynamics will be identical for all times to the original one. 

Next, we employ a chain of monotonicities, based on the graphical construction of Section 
2.3. Since we are after an upper bound on the set of minus spins, we can freeze to "— " all spins 
whose vertical coordinate is below L/i(4^). Therefore, we have just a dynamics in the set 

Y := [-4L^ + 1, 4Le - 1] x [L/i(4e), oo). 

In principle, its initial condition is such that the spin at site {xi,X2) G y is "— " if and only 
if X2 G [Lh{Ai) , LhQ{xi / L)\. The problem is however that the function x i— t- max(/i(4^), /io(a;)) 
is not concave, which prevents to apply directly Corollary 3.5 later. By monotonicity, we can 
modify such initial condition by adding extra "— " spins: we therefore stipulate that at time t = 
the spin at site (xi,X2) is "— " if and only if X2 G [Lh{Ai) ^ Lh{xi / L)]. Recall that h{x) > ho{x), 
so monotonicity goes in the correct direction. With some abuse of notation, we still call {a{t))t>o 
the dynamics thus modified and Aiit) the set of minus spins. We need a final step in order 
to map the evolution into the zero-range process. Note that, at time t = 0, the boundary of 
AL{t = 0), intersected with the strip [— 4L^ -|- 1/2, 4L^ — 1/2] x M, can be identified with the 
graph of a cadlag function 

Hl{-,0) : [-4L^ + 1/2, 4L^ - 1/2] ^ [Lh{4.C) - 1/2, oo) n Z, 

which is constant on intervals [n, n -|- 1) with n G Z and takes boundary values L/i(4^) — 1/2 at 
the two endpoints {Hl{x,0) is just a discretized version of Lh{x/L)). However, for time t > it 
is not true in general that the boundary of Aiit) is still the graph of a function, simply because 
the set Ahit) can be non-connected (see Figure 4). Let {(y^'^\t))t>Q be the dynamics obtained 
by erasing all the updates that would make AL{t) non-connected. It is easy to realize that, since 
Hl(-,0) has a single change of monotonicity (from non-decreasing to non-increasing, recall that 
h{-) is concave) such erased updates can only correspond to a "— " spin turning into a spin 
(see again Figure 4). Therefore, the set of minus spins of the dynamics {cr^'^^ {t))t>o dominates 
stochastically ^^(t): more precisely, we have shown that the coupling given by the graphical 
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construction implies that, w.h.p. and for all t > 0, 

ALit) C A^^\t) := U C (4.35) 

We let 

Hl{; t) : [-ALi + 1/2, 4Le - 1/2] ^ [Lh{A^) - 1/2, oo) n Z, 
denote the piecewise constant (random) function whose graph in the usual coordinates system 
(61,62) is the intersection between dA^l\t) and the strip [-4L^ + 1/2, 4L^ - 1/2] x M. Note 
that HL{-ALi + l/2,t) = HL{4L^-l/2,t) = L/i(4^) - 1/2. 

Equation (4.18) is proved if one shows that for any e small enough, w.h.p. 

yHl{x, sL^) < h{x/L, {a - 6)e) for every x G {-^L, ^L). (4.36) 

It is clear from Remark 3.3 that the function Hl{-, •) follows the dynamics described in Section 
3.2, with generator (3.3) (here we identify the function i?i(-,i) with an element of ^4^l-i/2i see 
(3.7)). According to Corollary 3.5 one has for arbitrarily small rj > 0, w.h.p, for all x S {—£,L, ^L) 

^Hl{x, eL^) < 4>{x/L, (1 + \\dJi\\oo)-^£) + V (4.37) 
where ||9a:^||oo = sup[_4^ 4g] \dxh{x)\ and 4'{x,L'^t) is the solution of 

dt4>{x,t) =ld^4>{x,t)^ 

(^(-4e,t) =m,t)=hm (4.38) 
(j){x,0) = h{x) for every x E [—4^,4^]. 

The equation (4.36) is thus proved if one has 

4){x, (1 + \\d.Ji\\oo)~'^e) < h{x, (a - 6)e) for every x e [-^,^]. (4.39) 

Note that by Lemma 4.3 (which is applicable because the second derivative of h{-) = /lo(') is 
uniformly continuous in (—2^,2^)) one has, uniformly on [— ^,^] 

^{x, (1 + \\dMoo)-^e) = 0(x, 0) + |(1 + WdMooyOl^ix, 0) + o(e) 

= ho{x) + ^(1 + WdMoorHdlhoiO) + r{x))e + o(e) (4.40) 

where r{x) tends to for 2; — t- 0. Finally, using (4.33), if ^ is chosen small enough so that both 
r(x) and ||9a;/i||oo sufficiently smaller than 5, 

^{x, (1 + WdMoor^^) < hoix) - (q - 6/A)e. (4.41) 

On the other hand one has 

h{x, {a - 6)e) > ho{x) -{a- 5/2)e, (4.42) 

which ends the proof of (4.18). □ 

4.4. Lower bound: proof of (4.4) and (4.2). The proofs follow the same ideas as those of 
Section 4.3: we need to control the dynamics for different portions of the interface separately 
(around the poles and away from them) using the scaling limit results provided by Theorems 
3.4 and 3.2. 
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4.4.1. Proof of (4.4). Equation (4.4) is absolutely crucial to start the proof of (4.2) and quite 
independent of the rest. The proof is very similar to that of [2, Theorem 2], so we only sketch 
the main steps. Set 

D ■.= {xe : d{x, (1 - 6)LS)) < 1} , D' := D ((1 + S'^)LS)y 

and consider a modified dynamics {a{t))t>o (whose law is denoted P), with the same initial 
condition as {a{t))t>o and the rules that: (i) after each update, any "— " spin which has more 
than two "+" neighbors is turned to "+", and the operation is repeated as long as such spins 
exist; (ii) the dynamics stops at the time f/) £>/, the first time when there is either a "+" spin 
in D or a "— " spin in D' . We define also td the first time when there is a "+" spin in D and 
Td,D'iTd the analogous random times for the original dynamics. 

Note that, by (4.3), w.h.p. td^di = td. Note also that the two dynamics can be coupled in a 
way that td^d' = td implies fu,D' = t^d 1^ td (thanks to point (i) above, since before ffj jji the 
modified dynamics has less "— " spins that the original one). Therefore, 

P(rD < £^2) = F{td < eL'^;TD,D' = td) + o(l) < P(fD < eL^; fz),D' = td) + o(l) 

and it suffices to prove for instance that 

i*{rD,D' < eL^^TD = fD,D') < exp(-7L) 

for some e = £{6) > 0,7 > 0. For this, one first observes (as in [2, Eq. (8.6)]) that when 
■^D,D' = T^D the difference between the number of "+" spins at time td and the number of "+" 
spins at time is at least cS'^Lp' deterministically, for some c > 0. 

Finally, (as in [2, Eq. (8.10)]) one proves that 

P{\{x : a^ieL") = +}| - \{x : a,(0) = +}| > cS^L") < exp(-7L) 

if e = e{5) is small enough. This is based on the fact (cf. [2, Lemma 8.5]) that, for times smaller 
than td^o', the rate of increase of the number of "+" spins is uniformly bounded by a constant. 

4.4.2. Scheme of the proof of (4.2). Given some fixed 5 > 0, we want to prove that for e > 
small enough, w.h.p. 

{1- {a + S)£)L& cAUeL'^), (4.43) 

or equivalently 

a^(eL^) = - for every x e {1 - {a + 5)e)L&. (4.44) 

Given ^ small enough (depending on 5) and u small enough (depending on ^), we define (cf. 
Figure 8) 

U ■.= {!- u)&, 

Ai(e) := [((1 - (a + 5)e)^) \ U] n [C, +00)^, (4.45) 
i?i(e) := [((1 - (a + 5)em \ U] n ([-^, i] x E+). 

and Ai, Bi, i = 2, 3,4 as the images of Ai{e), -Bi(e) by the rotation of angle i^. 
One has 

(l-(a + 5)e)^ = [/U (^U^*^ U (^U^*^ ' (4-46) 
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Figure 8. The large droplet is S! and (1 — (e(a + 5)))i^ is obtained by removing 
the external dark layer. The white central region [/, together with A\,B\ and 
its rotations (deformed rectangular regions) form a partition of (1 — eia + 5))^. 



and hence (using rotational symmetries), to prove (4.43), it is sufficient to prove that for e small 
enough, w.h.p. 

LUdAL^eL^). (4.47) 

LA^(e) C ^L(eL'), (4.48) 

LBx{e) C AL{eL^). (4.49) 

The first line, i.e. Equation (4.47), is a direct consequence of (4.4) provided that e is chosen 
small enough (how small depending on z/). Actually, one has the following stronger statement 
that will be useful for what follows: if e is small then w.h.p. 

L [/ C AiitL"^) for every t < e. (4.50) 

The main work is thus to prove (4.48) and (4.49). 

4.4.3. Proof of (4.49). This is similar to the proof of (4.18), except that monotonicities will be 
needed in the opposite direction. 

Let h : [—2^, 2^] i— t- M be a concave, twice differentiable, even function such that 

h{x) = hQ{x), VxG[-C,?], 
h{x)<h^{x), VxG [-2^,-0 U(e,2C] 

where /io(") was defined in Section 4.3.2 to be the graph of dSi n (M x ]R+) in the (ei,e2) 
coordinate system. Once ^ is fixed, we choose u and h such that the point (2,^, /i(2^)) lies in the 
interior of U . 

Using equation (4.50), we can freeze the spins with vertical coordinate Lh{2^) and horizontal 
coordinate in (— 2L^, 2L^) (we assume for notational convenience that 2L^ and L/i(2^) are in 
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Z*) to their initial value "— ", and w.h.p. , the dynamics we obtain is identical to the original 
one up to time eLp' . 




Figure 9. The sites in the dashed vertical lines are frozen to "+" and those of 
the horizontal bold segment to " — " so that the dynamics in the colored rectangle 
is independent of the rest of the system. At time i = the sites in the dark- 
colored region (whose upper boundary is determined by h{-) are "— " while those 
of the light-colored one are "+". The function /i(-) is such that the base of the 
dark-colored region is in LU . 



Next we use a chain of monotonicities based on the graphical construction of Section 2.3. 
Since we are after a lower hound on the set of minuses, we can freeze to "+" all the spins with 
horizontal coordinate ib2L^ and vertical coordinate larger than Lh{2^). Once this is done, we 
are reduced to considering the dynamics restricted to the set 

Y2 := [-2Le + 1, 2Li - 1] x [Lh{2i) + 1, oo), (4.52) 

as spins on its boundary are fixed. In principle, the initial condition one should consider is such 
that {xi,X2) E Y2 has spin "— " iff X2 G [Lh{2E,) + l,Lh^{xi/L)], but again by monotonicity, 
we can add extra "+" spins: we stipulate that, at time t = 0, {xi,X2) has spin "— " iff X2 G 
[Lh{2£^) + 1, Lh{xi/ L)]. With some abuse of notation, the dynamics thus modified is still called 

As for the proof of (4.18), we need a final step to map the dynamics onto the interface 
dynamics of Theorem 3.4, the problem being exactly the same as then: it is not true that 
the boundary of Aiit) stays connected for all t. The solution adopted in the previous section 
(leading to the dynamics {a^'^\t))t, see discussion before (4.35)) does not work here as we are 
now looking for a lower hound. 

Let {a^^\t))t be the dynamics that evolves like {a{t))t except that any spin that has three "+" 
neighbors is turned instantaneously to (see Figure 10). The coupling given by graphical 
construction implies that 

U C,=:Af{t)ciAL{t). (4.53) 

(3) 

Moreover our choice of initial condition guarantees that A^ {t) stays connected for all time, 
since the set ^ is convex. 

We denote by HL{-,t) the cadlag function [—2,^,2^] 1— t- M whose graph corresponds to the 

intersection between dA^^\t) and the vertical strip [— 2L^ + 1/2, 2L^ — 1/2] x IR- Note that 
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Figure 10. Light-colored (resp. dark-colored) squares denote "— " (resp. 
spins. In our modified dynamics a^^\ when a spin has three "-|-" neighbors, it 
is instantaneously turned to On the figure, if spin at A is updated and 

turns to "-|-", then the spin B has three "-|-" neighbors and therefore also turns 
instantaneously to 



Hii-jt) can be visualized as a collection of columns of width 1 and integer height. With this 
notation and (4.53), equation (4.49) is proved if one has w.h.p. 

yHl{x, L^e) > h{{a + 6)e, x/L)) for every x G (-^L, ^L). (4.54) 

Now we want to relate the dynamics of to that of Theorem 3.4. The relation is almost 
identical to that discussed in Remark 3.3, except for a slight difference in the way particles of 
types A and B annihilate in the zero-range process. Given Z 9 x = —2L^ -\- 1/2, . . . , 2L^ — 1/2, 
we say again that there are n > particles of type A at time t at site x if limj^_^^+ Hi{y,t) — 
liniy^x- Hiixjt) = n and that there are n > particles of type B if the same difference 
equals —n. Then it is easy to realize that, under the dynamics {a^^\t))t^Q, each particle 
performs a symmetric simple random walk with jump rate l/(2?i) both to right or left (with 
n the occupation number of the site where the particle is), and that particles of different type 
annihilate immediately if they are at sites of distance 1 (and not on the same site): this is 
the effect of flipping instantaneously "— " spins with more than two "+" neighbors. Note also 
that, due to convexity of /i(-), particles of type A are always to the left of particles of type B. 
Therefore, if we take i?L(-, t) and we eliminate one of the columns of maximal height (see Figure 
11) (note that there are always at least two), the modified height function thus obtained follows 
exactly the evolution of Theorem 3.4. 




N ~ 1 



Figure 11. Left: the height function associated to the "+/— " boundary for 
the dynamics a^^\t). Right: the same height function, with one of the highest 
columns removed; this follows the same evolution as in Theorem 3.4. The fact 
that the new interface is step shorter makes no difference in the macroscopic 
limit. 
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Of course, the erased column does not change the scahng hmit so that one can apply Theorem 
3.4 and Corollary 3.5 and get that for any t and r] > 0, w.h.p. 

^-Hl{x, L^t) > ^{x/L, t)-r] for every x e {-2^L, 2^1} 



L 



(4.55) 



where 

dt4>{x,t) = ^d'^(p{x,t) 

m,t) =^i-2C,t) = hi20 (4.56) 
(p{x,0) = h{x) for every x G [—2^,2^]. 
Therefore, (4.54) is proved if one can check that 

^(x, e) > /i(x, (a + (5)e) for every xG[— ^,^]. (4-57) 

The above equation is proved is the same manner as (4.39): one just needs to choose ^ small 
enough. 

□ 

4.4.4. Proof of (4.48) . First of all, one freezes to " — " all the spins on the cross-shaped region 
of sites in LU (cf. (4.45)) such that at least one of their coordinates is ±1/2. Equation (4.50) 
guarantees that if e is chosen small enough, w.h.p. the so-obtained dynamics coincides with the 
original one up to time eL^ if e is small enough. 

Then one defines {a^^Ht))t ^ as the dynamics obtained by changing the initial condition in 
the following manner: all spins {x,y) S with either |a;| ^ L{1 — v) or \y\ ^ L(l — v) are 
changed from "— " to "+" (recall that v is the constant that enters the definition (4.45) of U) 
and therefore they stay "+" forever, since they have at least three "+" neighbors. Note that, 
this way, the evolution in each quadrant of (Z*)^ is independent. By monotonicity, we get that 
w.h.p, for every t < eL^, 

U C,=:A^^\t)cALit). (4.58) 

and therefore (4.48) is proved if one can show that 

LAi{e) cAf{eL^). 



(4.59) 



Next, note that dA!"^\t) n in the coordinate system (fi,f2) is the graph of a random 
piecewise linear function 

[l-u)L {l-u)L 



Fl ■■ 



V2 V2 

which undergoes the corner-flip dynamics described in Theorem 3.2 (apart from space rescaling 
by a factor \/2)- For this reason on gets that, w.h.p.. 



lim sup sup 



jFL{xL,tL^)- g{x,t) 







where 



dt9{x, t) 



ld^g{x,t) 



g{x, 0) = f{x) for every x G 

and / is the profile of the initial condition, i.e. 

fix) :=mm{fo{x),{l-u)V2-\x\)) 



^/2 ' V2 



(4.60) 



(4.61) 



(4.62) 
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Let p (resp. pi) be the point on dS> whose coordinates {x^y) (resp. (xi,X2)) in the coordinate 
system (ei,e2) satisfy x > = 1 — v (resp. xi > 0, yi = /io(0)- Call —d < (resp. —di < 0) 
the horizontal coordinate of p (resp. of pi) in the coordinate system (fi, £2) . 
In view of (4.60) and of Definition (4.45) of Ai{e), equation (4.59) is satisfied if 

g{x,e) > f{x,{a + 5)e) for every x£{—di,di), (4.63) 

The proof of this is very similar to that of (4.24) provided that that / coincides with /q in 
a domain containing strictly {—di,di) (this guarantees for instance that 9^/(-) is uniformly 
continuous in a domain containing (— c?i,di), so that the drift dtg is continuous in time, cf. 
Lemma 4.3). For this to hold, it is enough to assume that d > di, i.e. that i' in (4.45) has been 
chosen sufficiently small as a function of ^ so that 1 — > /lo(^). □ 

5. Proof of Theorem 2.1: existence of anisotropic curve-shortening flow with 

CONVEX initial CONDITION 

Let us first recall some properties of the support function h{-) of a convex curve 7. First of 
all, if 7 C 7' then h{6) < h'{9) for every 6. Next, the support function is related to the curvature 
and to the length -^(7) of 7 by (cf. [12, Lemma 1.1]) 

dlh{e) + h{e) = ^^ (5.1) 

2n r2n ^ 



U;) = j^ hm» = j^ ^dO. (5.2) 

Also (cf. Lemma 4.1.1 in [13], with the warning that what they call \s 9 — 7r/2 for us), the 
Cartesian coordinates {x{6),y{6)) of the point of 7 where the outward directed normal forms an 
anticlockwise angle 9 with the positive horizontal axis can be expressed as 

x{9) = h{0) - f ^^ds (5.3) 
Jo Hs) 

y{9) = hin/2) + [' ^ds. (5.4) 

Under the flow (2.7), the time derivatives of area and length are (cf. [12, Lemma 2.1]) 

— Area(7(t)) = - / a{9)d9 (5.5) 
dt Jo 

L(7(t)) = - / a{9)k{9,t)d9. (5.6) 







dt 

For the moment these are formal statements since we do not know yet that the flow exists. 



5.1. Proof of Theorem 2.1. Uniqueness of the flow is trivial, so we concentrate on existence. 
First of all, we need to regularize the functions a(-) and k{-). Given < w < 1 we define a^'"\-) 
to be a family of smooth approximations of the anisotropy function a(-). More precisely: 

Assumption 5.1. 

(1) a("')(-) is liT-periodic and C°°; 

(2) a^'^\9) a{9) uniformly in 9; 

(3) for fixed 9, the function w 1— >• a^^\9) is non-increasing; 

(4) the function a^^\-) is Lipschitz, uniformly in w > (this is possible because the function 
a(-) itself is 1 — Lipschitz); 
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(5) the functions w i— )• HSgO^'^) ||oo := max0\dga^'^\9)\ and w i— ||5|a("') ||oo are bounded, 
uniformly for w in any compact subset of (0, 1). 

A possible choice is 

where g^""^ is a centered Gaussian of variance w . In the convolution it is understood that a(-) 
is seen as a 27r-periodic function on M and is chosen so that a^'^\-) satisfies the monotonicity 
with respect to w. It is easy to check that one can choose = —Cw for some suitably large C. 
Indeed, monotonicity in w is guaranteed if for w' < w one has 

(9^") -9^"''))||oo. 

On the other hand, since o(-) is Lipschitz, one sees easily that ||a*((7("'-' — 9^"''-')||oo = 0{w — w'). 

Also, we approximate 7 with a sequence of convex curves {'J^^^)o<w<i that satisfy the following 
properties: 

Assumption 5.2. 

(1) 7^"') D 7('^') D 7 or equivalently h^'"\-) > > h{-) if < w' < w; 

(2) lim^_>>o = h{-) uniformly in 9, so that 7 is the limit 0/7^"'^ in the topology of the 
Hausdorff distance; 

(3) the Lipschitz constant L(A;("')) of the curvature function k^'^\-) is finite uniformly in w, 
k^'^\-) — ;> k{-) uniformly and limsup^^QL(A;('")) < L(A;); 

(4) The three first derivatives with respect to 6 of k^''"\9) are bounded uniformly for w in 
any compact subset of (0, 1). 

(Like for the regularization of a(-) into a^'^\-), a possible construction of h^'^\-) is obtained 
convolving /i(-) with a Gaussian of variance w'^ and adding a suitable constant e^). 
For the regularized mean curvature motion, it follows from [12] that the equation 

idth^^\9,t) = -a^^\9)k^^\9,t) 

|M"')(e,0) = h^'"){9) ^ ■ ' 

admits a solution corresponding to a flow of curves {'^^'^\t))t>Q which remain convex and shrink 
to a point in a finite time 

if ■= ty' = Area(7("')(0))/ J a^'"\9)d9 

(cf. (5.5) with a(-) replaced by a^^\-)). For lightness of notation, we will often write h{-, •), 7(t), a(-), 
etc. for the regularized quantities h^^\- , ■) , ^^^\t) , a^'^\-) , etc. Also, we abusively identify 
a convex curve 7 with the set enclosed by it. Thanks to Assumption 5.1, we have that 
f^^a^'"\9)d9 /o^^o(0)d0 = 2 as -uj ^ and therefore t'-f^ = tf{l + o(l)) with tj defined 
in Theorem 2.1 when w ^ 0. 

From (5.1) and (5.7) one can check that the curvature satisfies the parabolic equation 

(dtk = k''dl{dk) + ak^ 

\k{0,0) = k{9). ^ ■ ' 

Also, following [13] it is possible to see that the curvature function stays C°° until tj (since a is 
C°°). However, estimates on the regularity will not be necessarily uniform in the regularization 
parameter w and we will need to be very careful on this point. 
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For fixed t, let 

7(t)= lim7(-)(t) (5.9) 

where convergence is in the Hausdorff metric. A posteriori, since we will see that {^{t))t provides 
the (unique) solution to our curve-shortening equation, it follows that the limit (5.9) does not 
depend on the choice of regularization. Existence of the limit is guaranteed by the fact that 
7^^" ^ (t) C 7^^"^ (t) if w' < w and t < t^^ (this follows from the fact that a^"') {0) is decreasing 
in w and that the curve is smooth at all times). One has to use the Blaschke selection theorem 
[8, Th. 32] which says that a family of convex subsets of a bounded subset of admits a 
sub-sequence converging to a non-empty convex set. Convergence in Hausdorff distance also 
holds for the boundary curves. 

Since the volume is continuous in the topology induced by the Hausdorff metric [8, Ch. 4] 
we also see that Area{-f{t)) = Area{-f) - t a(6')d6l = Area{j) - 2t; for t tf the curve 7(i) 
shrinks to a point (its diameter shrinks to zero). 

Definition 5.1. Fort < if /ei femax(i) (resp. femin(i)J be the maximal (resp. minimal) curvature 
of ^{t). We let fcmax := A;max(0) and similarly for fcmin and amax(min) max5i(mine)a(6'). Also, 
^min(max) 0"-^? amin(max) o,re defined similarly to k^hiimux) , "'mm{max) but with A;(-),a(-) replaced by 
k{.),a{-). 

It is crucial that kmaxit) stays bounded, uniformly for w small, as long as the disappearance 
time is not approached: 

Proposition 5.2 (Regularity estimate). Assume that the curvature function k{-) is Lipschitz. 
There exists wq > such that, for every b > 0, t < tf{l — b), < w < wq one has 

^max(t) < Ci (5.10) 

and 

max\dg{d{e)k{e,t))\ <C72 + L(A;) (5.11) 

6 

where we recall that L(A;) is the Lipschitz constant of the function k(-). The constants Ci and 
C2 depend only on b and on fcmax- 

We will prove 

Theorem 5.3. The flow of curves (7(t))j<j^ defined in (5.9) is a classical solution of the 
anisotropic curve shortening flow (2.7) for < t < tf . 

Proof of Proposition 5.2. The proof is based on ideas of [12]. However, it is important to make 
sure that estimates are uniform in w < wq (in [12] the anisotropy function a(-) is assumed to be 
C^, so there was no need to regularize it). 

Fix w > {). First we get a lower bound on A;min(i)- Note first of all that at time zero the 
minimal curvature is bounded away from zero (uniformly in w): indeed, using (5.2) and the fact 
that the curvature function is L(A;)-Lipschitz, 

^(7(0)) = r T^d. > 2 r , \,..A0 = ^ log H^t±^. (5.12) 



JO k{9) - Jo K^in + Uk)0 Uk) " k 

Then, since the length of 7(0) is finite, k^^[^ must be positive. 
Set for simplicity 

g = g{e,t) = a{e)k{0,t). 
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Formula (5.8) gives 



dtg = \{g^dlg + g^) =: g{e,t)u{e,t). (5.13) 



This, together with the fact that a(-) and k{-,t) are smooth, imphes that 



d mine/ . . 

— min5 > ^ > (5.14) 

at flmax 

(at the minimum point the second derivative is positive) so that 

^min(i) ^ ^min ^ C'^min > (5.15) 

'^max 

with C independent of w (say for w < wq) thanks to the uniform convergence a^'^\-) — t- a(-) 
and A:("')(-) ^ A;(-). 

Next the real work: bounding kmaxit) uniformly in w. From (5.13) one sees that, since a(-) 
and k{-,t) are smooth, 

— max (7 < x (max (7)^. (5.16) 

From this one immediately gets that /Jmax 

(t) is upper bounded uniformly in < wq, up to 
some time ti depending only on fcmax- However the solution of x = explodes in finite time, 
certainly before the time tf when the curve shrinks to a point, so we need to do better. 

For this, we define z{t) = mi'n.Qu{9,t) (cf. (5.13)). Then, taking the derivative of u with 
respect to t shows (cf. Lemma 4.2 of [12] for details) that 

^z(t) > 2z{tf 

so that if z(0) > we get z{t) > 0, if z(0) < we get z{t) > -l(l/|z(0)| + 2t). Altogether, we 
get that 

n(0,t)>-l (5.17) 

uniformly in 9 and w < wq. Now we use this to get a uniform bound on 1 1 ^0^1 1 00 tGrnis of 
^max(^)- Without loss of generality suppose that there exists Bi such that dQg{9\^t^ — ||f^(?5'l|c)o 
(if this is not the case one can still find 9i such that dgg{9i,t) = —\\dgg\\oo and apply the same 
method). Let also 62 > Oi be such that dgg{62,t) = (such angle exists since g is periodic). 
Then, from the definition (5.13) of u, 



|9,9||oo = - r digde = - r f _ a(0, t)k{e, t)] dd 

J 01 J 01 \k[d,t) J 



.k{d,t) 

^ ^ + - 0lXaxfemax(i) < ^^^^^ + Cji,...{t). (5.18) 

^tJoi k{6,t) t 
In the last inequality we used (5.2) and then (5.6) which says that L{'y{t)) < L{'y{0)){< 
2L(7(0))). Since g = ak and by assumption o is and Lipschitz uniformly in w, one de- 
duces that 

\\dek\\oo < ^^^^ + C5^w(t) < C6(t)fcmax(t) (5.19) 
and Ce can be chosen to be decreasing in t. From this it is trivial to see that, if 9o is such that 
k{eo,t) = femax(i); one has 

k{0, t) > Imax(t)/2 whenever 16* - 6*0! < a{t) (5.20) 
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for some a{t) increasing in t (it could vanish for f — )• 0). Next, one proves that for t < (1 — b)tf 

Eit) := r d{9) logigie, t))de < Cj (5.21) 
Jo 

where C7 depends only on Omax and on b and on the maximal curvature fcmax of the initial curve 
7(0). Indeed, (5.21) is obvious for t = 0, since the initial curvature is bounded by assumption. 
To get the control for t > 0, one observes (cf. Propositions 5.3 and 5.4 of [12]) that 

The prefactor is bounded since b > and the time-integral of the time-derivative of the length 
gives at most ^(7(0)). At this point we are almost done: using (5.20) 

Cj> r ~a{9)^ogig{e,t))de (5.22) 
Jo 

> 2a(t)dmin log(amin^max(i)/2) + 27ramax log[min(l, a^inkminit))] (5.23) 

and this (recall that k„im{i) > C'^min > 0, cf. (5.14)) gives us an upper bound on A:max(i) 
uniformly in w < wq and t < (1 — b)tf: up to ti one uses the upper bound which comes 
from (5.16) and after ti the one from (5.22); Eq. (5.10) is proven. When t approaches the 
disappearance time if (i.e. when b approaches zero) the upper bound diverges (because C7 
diverges), as it should. 

Equation (5.19) says that the curvature function is Lipschitz with a Lipschitz constant C that 
depends on t,b and L(0) but not on w. This is not yet the desired (5.11) because the bound 
diverges for t — t- 0. To prove (5.11) remark that, using (5.13), 

dtdea = de (Cdh +C]=-^ {g'dh + (5.24) 

+i {2gdegdlg + g^dlg + 3g^deg) . (5.25) 
a 

At the point where d^g is maximized, dgg cancels and d^g is non-positive. This, together with 
the boundedness of g uniformly in w < wq, 6 S [0, 27r] and t < {1 — b)tf, implies 

dtmsixdgg{e,t) < Cs{l + max dgg {6, t)). (5.26) 

6 

where Cs just depends on A:inax and b. Integrating with respect to time, one gets 

maxd0g{9,t) < m&xdeia^'"He)k^'"\e)) + Cg 

6 
with C9 depending only on Cg. Also, observe that 

dg{a^^\e)k^^\e)) < {?,/A)\dgk^^\e)\ + Cio < (3/4)L(A:("')) + Clo 

with Cio a constant depending on /cmaxj since for w small a^lx < (3/4) and a^'") is uniformly 
Lipschitz. Finally, from Assumption 5.2 (3), we can conclude de{a^'"\9)k^'"\e)) < Cio + L(A;) 
for w small. An analogous lower bound can be found on dtming dgg{6,t) and this gives (5.11). 

□ 

Following [13] it is possible to prove that, once we have bounds on curvature and on ||t^^?^(*, ^) ||oo ; 
for every n > 2 and t < tf{l — b) the derivatives dgg{6,t) are also bounded. The bounds we get 
are in general not uniform in w but this is not very important for our purposes. Indeed, we will 
need only: 
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Proposition 5.4. Fix b > 0. There exists a function c{w), which is non-increasing with respect 
to w £ (0, wq] such that for t < {1 — b)tf 

max\df h{e, t)\ <c{w). (5.27) 

6 

Proof Recall (5.7) and (5.13): 

drh = -l{g'dlg + g'). (5.28) 
a 

Thus we just have to bound dgg, since we have already proved that g itself is bounded. For this, 
we adapt the method used by Gage and Hamilton in [13] for the special case of the isotropic 
curve shortening flow where a = I. What they observed [13, Lemma 4.4.2] is that, if the 
curvature and its ^-derivative are bounded (which we proved in Proposition 5.2), the t-derivative 
of ^{t) := J^^[dgg{6,t)]'^d9 can be upper bounded by a constant times ^{t) itself and then one 
can integrate the inequality with respect to t to get a bound on <I>(t) in terms of <I>(0). In 
our case, with a similar computation, we find that {d/ dt)^{t) is upper bounded by ^{t) times 
a constant depending on H^eo^"'^ ||oo, which is finite uniformly for w < 1. Since ^*(0) is also 
bounded for w in compact subsets of (0, 1) (cf. Assumption 5.1 (5) and Assumption5.2 (4)), we 
get that <I>(t) < ci{w) for w G (0, 1) and t < (1 — b)tf and we can choose ci to be decreasing. In 
general, ci will diverges when w approaches zero. 

A similar computation (cf. [13, Lemma 4.4.3] when a{0) = 1) shows that 

^(0 := r[dlg{e,t)fde<C2iw) 
Jo 

with C2(-) decreasing in w £ (0, 1). Then one uses the fact that for a smooth, 27r-periodic 
function / one has (cf. [13, Corollary 4.4.4]) 

\\f\\l<c r{f + {defr)de 

Jo 

for some universal constant C, applied with /(•) = dgg{-,t), to get that II^^^Hoo < cs{w) as we 
wished. 

□ 

Proof of Theorem 5.3. We are now ready to prove that (7(t))t provides a classical solution of 
(2.7). This is based on the following easy consequence of the Arzela-Ascoli Theorem: 

Lemma 5.5. Let /("^ be a sequence of periodic functions on [0, 27r], such that both sequences 
Z*-"^ and dxf^^^ are uniformly bounded and equicontinuous. If f^^^ f as n ^ oo, then f is 
and dxf = liinndxf^^\ where the convergence is uniform and does not require sub-sequences. 

First of all, we note that h{-,t) does converge (for w ^ 0) to h{-,t) for every fixed t < tf. This 
just follows from the fact that 7(t) converges to 7(i) in terms of Hausdorff distance. Applying 
Lemma 5.5 and recalling (5.1), we get that, for t fixed, deh{6,t) and k{6,t) converge to deh{6,t) 
and k{t,9) respectively and that convergences are uniform in 9 (knowing that the curvature 
is Lipschitz is important here). Note by the way that k{-,t) is Lipschitz, since ||c)eA;(-, t)||oo is 
uniformly bounded. 

Then applying dominated convergence (which is allowed in view of Proposition 5.2), one gets 
that 

h{e,t)-h{9,s) = - f a{9)k{9,u)du, (5.29) 

J s 
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which is an integrated version of (2.7). To get the stronger statement (2.7), we need to prove 
that k{6,t) is continuous function of t. 

First of all, we prove that one can find a function e : (0,1) 3 w ^(^) ^ increasing and 
going to zero as w — )■ such that for all 9, for all t < (1 — 6)t/, 

\k{e,t)-k{e,t)\<e{w). (5.30) 

If this were not the case then, thanks to the fact that k{-,t) and k{-,t) are uniformly Lipschitz, 
we would have, say, for arbitrarily small w and for some e > 0, 

k{t,e) - k{t,9) > e 

for ee[e,e + e\ for some 6 G [0, 27r]. But then, since (cf. (5.1)) 

1 1 



{d'e + i){h{e,t)-h{e,t)) 



this would contradict the uniform convergence of h{-,t) to h{-,t). 

On the other hand, from Proposition 5.4, for all 9 and for all t, s < {1 — b)tf 

\k{e,t) -k{e,s)\ <c{w)\t-s\. (5.31) 
Together with (5.30) this implies that 

\k{9, t) - k{6, s)\ < inf {2e{w) + c{w)\t - s\) . (5.32) 

The right-hand side clearly tends to zero with |t — s| (choose a sequence {wk} tending to zero. If 
c{wk) does not diverge we are done. Otherwise, compute the right-hand side for the w = Wk with 
the largest value of k such that c{wk) < \t — s|~^/^). This shows that t i— )■ k{9,t) is continuous 
away from tj and the proof is complete. □ 

6. Proof of Theorem 2.2: evolution of a convex droplet 

The proof is very similar to that of Theorem 2.3 in the scale-invariant case (Section 4), and 
therefore it will be only sketched. We will also try to use as much as possible the same notations 
as in Section 4. 

First we present two statements that are analogous to Propositions 4.2 and 4.1: 

Proposition 6.1. Let T> he convex with a Lipschitz curvature function. For every a > 0, w.h.p. 

AiiL^t) C LV^'^^ for every t>{) (6.1) 

(recall definition (2.8) ). Moreover, for every a > there exists ei(a, fcmax) > such that w.h.p 

AULh) D for every te [0,ei]. (6.2) 

Proof. The proof of (6.1) is essentially identical to that of (4.3), so we give no detail. As for 
(6.2), given a it is possible to give a finite collection such that: 

• each is an open convex subset of M^, obtained from (the interior of) the invariant 
shape & via a suitable translation and shrinking; 

• ^^i C V for every i; 
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Given rj > 0, thanks to Proposition 4.2 there exists e > such that, w.h.p., for every t < e one 
has 



Here we use monotonicity (because &i C T>) and the fact that the union of a finite number of 
events which occur w.h.p. still has probability tending to 1. Note that the choice of e is depending 
on f] but also on the diameter of the smallest set in the collection and consequently on 

^max- Then, if rj is small enough (depending on a) it is clear that Uji^^ ''^ D (recall that 

the are open sets, so that every x E D(^") is contained in the interior of at least one □ 

Proposition 6.2. For all 6 > there exists eo(<^) ^mim ^max) > such that for all < e < Eq, 
w.h.p., 

AL{L^e) C LV{e{l - 6)), (6.3) 

and 

ALiL'^e) D LV{e{l + S)) (6.4) 
where we recall that T)[t) is the set enclosed by the curve 7(t). 

Proof of Theorem 2.2 assuming Propositions 6.2 and 6.1. It is enough to prove (2.9) for t < 
(1 — b)tf for arbitrary 5 > 0. Then, the statement for t > (1 — b)tf and also (2.11) follows from 
the fact that the disappearence time of a droplet of diameter i is w.h.p. 0{£'^) (recall that 7(i) 
shrinks to a point t — t- in the sense that its diameter converges to zero). Define fcj^jj^ > (resp. 
^max < to be the infimum (resp. maximum) of A;min('S) (resp. /cmax('S)) on [0, (1 — b)tf]. Fix 
6' small and let e < eo('5'; ^min' ^max) ^ < ^li^f^i ^nia.x) with eo,ei defined in Propositions 
6.1 and 6.2. Using the Markov property and the monotonicity of our process we get that, w.h.p., 
for any k such that ek < (1 — b)tf 

AUL^ke) C LV{ke{l - 6')). (6.5) 
From (6.5) and Proposition 6.1 we get that w.h.p., for every t ^ (1 — b)tf, 

-\ e{l - 5')J J C L [D{{t - e)(l - 6'))] ^"^^ . (6.6) 

Setting e' = tf6' + e this implies that w.h.p. 

AUL^t) C LV{t - e')^^/^^ for every t ^ (1 - b)tf. (6.7) 

Finally observe (this follows from (2.7)) that the Hausdorff distance between 'D{t — e') and I^(t) 
is at most e'k*^^^maxQ \a(9)\ so that if e' is chosen such that 

e'k*^,,max\ai9)\<6/2. (6.8) 

u 

we get (2.9). 

The lower bound is proven similarly and this is where one has to use the assumption e < 

ei{6/2,k*^,J. □ 

6.1. Upper bound: Proof of (6.3). 

Definition 6.3. Define {Pi{t))j^^ to be the four "poles" of'D{t), where the tangent vector is 
either horizontal or vertical (recall that T>{t) is strictly convex at all times under our assump- 
tions, cf. discussion after (5.12), so that the four poles are distinct and uniquely defined). 
Pi{t) denotes the "north pole" and the others are numbered in the clockwise order. Denote by 
{x{Pi{t)),y{Pi(t))) (resp. {u{Pi{t)),v{Pi{t)))) the coordinates of Pi(t) in the coordinate system 
(fi,f2) (resp. (ei,e2)j. When t = we omit the time argument. 



Al{LH) C L 



V 
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Figure 12. The larger convex set is T> and the smaher one is — (5)). The 

poles Pi of T) are marked with black dots (for convenience we have chosen Pi one 
the vertical axis and P2 on the horizontal one). The graph in (fi, £2) of the anti- 
clockwise portion of dT> between A and B is /(-jO) and the graph in (ei,e2) of 
the portion of dT> between P4 and P2 is 0). For the proof of (6.11), boundary 
spins to the left of li are set to "— " below Pi and "+" above; boundary spins 
below l2 are set to "-" to the left of P2 and "+" to the right. 



An equivalent formulation of (6.3) is: for all 5 > and e small enough w.h.p. 

a^{£L^) = + for every x G L [P(e(l - 5))f . (6.9) 

Given some small ^ we divide ['C(e(l — 5)yf in eight pieces {Mi)^^^ and {Ni)f^-^ as follows (this 
is analogous to the definition (4.16) in the scale-invariant case, cf. Figure 6): 

Mi{e,0 := (KPi) + e,oo) x KP2) + 00)) \ P(e(l - 5)) (6.10) 

while A^i(e, ^) is the infinite component of {[u{Pi) — S^,u{Pi)+^]xM)\D{e{l — 5)) which contains 
Pi. The sets Mi,Ni are defined analogously for i = 2, 3, 4, so that [P(e(l - 5))]'' = [jt=iiMiUNi). 
Equation (4.15) is proved if one can prove that for every i, and e small enough, w.h.p. 

ax{eL^) = + for every x G LMi{e, ^) (6.11) 
o-^(eL2) ^ _^ fQj. g^gj.y ^ g LNi{e,£,). (6.12) 

Of course one can focus on i = 1, the other cases being obtained by a permutation of coordinates. 

6.1.1. Proof of (6.11). We use the notation f{-,t) for the function whose graph in the coordinate 
system (fi,f2) is the portion of &D[t) which goes in the anti-clockwise direction from point A 
where the tangent forms an angle 7r/4 with the horizontal axis (cf. Figure 12) to the point B 
where the angle is (5/4)7r. The domain of definition of f{-,t) decreases with time (because T>{t) 
shrinks) but for t small enough it includes [x(Pi), x(P2)]. Let Di be the "triangular-shaped" 
region delimited by dT), by the vertical line ii passing through Pi and by the horizontal line £2 
passing through P2 (note that Di may not be included in D). 
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We consider a modified dynamics in the north-east quadrant [Lu{Pi),oo) x [Lv(P2),oo) de- 
hmited by the hnes L£i,Li2. All the spins are initially "— " in LDi and "+" otherwise. As for 
boundary spins, the spins at distance at most 1 to the left of Lii are frozen to "— " if they are 
below LPi and to "+" if they are above. The spins at distance at most 1 below Li2 are frozen 
to "— " if they are to the left of LP2 and to "+" otherwise, see Figure 12. In the quadrant under 
consideration, this dynamics dominates the original one (for the inclusion order of the set of 
"— " spins). Let Fi{-,t) denote the function whose graph in (fi,f2) is the interface between "— " 
and "+" spins for this dynamics. Using exactly the same argument as in (4.22) we get that 



lim sup sup 

xelx{Pi),xiP2)] t ^ T 



jFL{xL,tL^) - g{x,t) 



(6.13) 



where g is the solution for t^O and x € {x{Pi), x{P2)) of 

dtg{x,t) = \d'^g{x,t) 
< g{;t) = f{;0) (6.14) 
^g{x{Pi),t) = y{Pi) and ^(^(^'2),*) = ^(^2). 

We are thus reduced to prove that for every xi,X2 satisfying x{Pi) < xi < X2 < x{P2) and 
every x € (xi, X2) 

g{x,e)<f{x,{l-5)e). (6.15) 

Lemma 4.3 (which is valid also in this case) allows us to write that for any fixed r/, for e small 
enough, 

g{x, e) ^ /(x, 0) + I {dlf{x, 0) + 77) . (6.16) 

We are left to estimate the right-hand side of (6.15). For any 9 G (0, 7r/2) and s > define 
x{9, s) to be the fi coordinate, in the (fi, £2) coordinate system, of the point of 7(5) where the 
outward normal vector forms an anticlockwise angle 6 with the horizontal vector ei. Note that 
for s ^ x{-,s) defines a bijective function. We denote 9{-,s) its inverse. 

It is more practical for the purposes of this section to rewrite the curve-shortening flow in 
the (fi,f2) coordinate system. Using the explicit expression (2.4) of a(0), some trigonometry 
and the expression |/"(rE)|/(l + (/'(a;))^)^/^ for the absolute value of the curvature at the point 
(x, /(x)) of the curve given by the graph of a function x 1— )■ f{x), one gets that for 9 G (0, 7r/2) 

a{9)k{9, s) = -\dlf{x{9, s),s) cos{9 - j) (6.17) 



and 



a{9{x,s))k{9{x,s),s) 1 3.. x 

= - cos(g(x,.)-f) = l9Jix,s), (6.18) 

so that 

f{x,{l-6)e) = f{x,0) + -dlf{x,s)ds. (6.19) 

We need therefore to prove time-regularity of 5^/(-, s): 
Lemma 6.4. For every x{Pi) < xi ^ X2 < x{P2) 

limsup{|(9t/(x,s) - (9t/(x,0)|,s G [0,t] and x e [Si,i2]} = 0. (6.20) 

t-~>-0 



36 



H. LACOIN, F. SIMENHAUS, AND F. L. TONINELLI 



Proof of Lemma 6.4- Recall from Section 5 that the curvature function k(9, s) is jointly contin- 
uous in {6,s). Thus using equation (6.18) it is sufficient to prove that 6{x,s) is a continous 
function in s uniformly in x: 

limsup{|6'(x,s) - 61(2;, 0)1, s G [0,t] and x G [5i,52]} = 0. (6.21) 
This comes from the continuity of x(0, •) 

limsup{|x(6i,s) - x{e,0)\,s G [0,t] and 9 G [0, ^]} = (6.22) 
and from the fact that x(-, s) is strictly monotone: for t ^ 0, 

mf{\dexie,s)ls^t,9e[0,^]}>0. (6.23) 
Both properties are a consequence of 

,n ^ / N r C0S(9' -TT/A)d9' 

x{9, t) = x{7r/A, t) - / ^ ' ' (6.24) 

Jn/i K[t) ,t) 

which is easily derived from (5.3)-(5.4). □ 
We finally get that for x G (xi, X2) and e small enough, 

fix, (1 - 6)e) > fix, 0) + (1 - 6)^idlfix, 0) - r?). (6.25) 

Thus, combining this with (6.16), (6.15) is proved if one has 

d'jix, 0) + ^ (1 - 5) idlfix, 0) - v) (6.26) 

i.e. 

2?7 + 59|/(x,0) ^ 0. (6.27) 

For this it is sufficient to have r] small enough, since (cf. (6.18)) sup{5^/(x, 0), x G [x(Pi), x(P2)]} 
can be upper bounded by a negative constant times the minimal curvature /cmim which is strictly 
positive. □ 

6.1.2. Proof of (6.12). Set /i(-,t) to be the continuous concave function whose graph in the 
(ei, 62) coordinate system is the portion of of 7(i) which goes from i-2(i) to -P4(t) with the anti- 
clockwise orientation. Given a small r] choose ^ small enough so that sup{\dxhix,0)\,uiPi) — 
e ^ X ^ uiPi) + C}<v- 

Consider /i(-) the function equal to /i(-,0) on [uiPi) — 2,^,n(Pi) -I- 2^] and linear outside. 
Assume for definiteness that /i(t((Pi) — 4^) < /i(ti(Pi) -1- 4^). Define ^~ = n(Pi) — 4^ and 
= inf{x > u(Pi),/i(x) = /i(^~)}. We consider the restriction of h to , C^] and still call it 
h. Define 

:= [e+,00) X [/i(^+),oo), J2 := (-oo,r] X [M^+),oo). (6.28) 

We consider the same chain of monotonicities as in the scale-invariant case (Section 4.3.2) and 
we end up with a dynamics in the half-strip [L(^~ , L^'^] x [L^(^"'"), 00) with boundary spins frozen 
to "4-" in L(J^ U J2)_and to "-" in Z* x (-00, L/i(^+)] and an initial condition with "-" spins 
under the graph of L/i(-/L). Also, the dynamics thus obtained does not allow moves that make 
the interface non-connected. Calling (cr2(i))oo this dynamics, (4.35) is satisfied. 

Define Hi : [L^~ ,L^~^] — >■ Z to be the function whose graph in (ei, 62) is the interface between 
"+" and "— " spins. We have to prove 

^-HLiLx,£L^) < /i(x, (1 - (5)e) for every X G (n(Pi) - ^, n(Pi) + ^). (6.29) 



L 
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Following the same steps as in (4.37) to (4.40) one finds that the left-hand side of (6.29) is upper 
bounded w.h.p. by 

h{x, 0) + |(1 + r?)-2(aXn(Pi), 0) + r{x)) + o{e), (6.30) 

where r{x) tends to when x — )■ u{Pi). 

To estimate the r.h.s of (6.29), one remarks that, in analogy with (6.18), 

dsh{x, s) = —a{6{x, s))k{9{x, s),s)/ sm{9{x, s)) (6.31) 

so that dth{x,t) is continuous in x and t (since 9 is around 7r/2, sin(6'(x,s)) is bounded away 
from zero). Moreover 

dth{u{Pr),Q) = ^5>(n(Pi),0), (6.32) 

which can be obtained directly from a(0) = 1/2 and from the fact that the curvature of P at 
the north pole Pi equals minus the second derivative of /i(x,0) computed at x = u{Pi). Thus 

h{x, (1 - 5)e) ^ h{x, 0) + (1 - (5)|(1 + r])dlh{ui,0), (6.33) 

and (6.29) is proven (combining (6.30) and (6.33)) choosing r] and ^ small enough. 

6.2. Lower bound: Proof of (6.4). We are confident that the reader is by now convinced 
that the proof of Theorem 2.2 is essentially identical to that in the scale-invariant case, modulo 
the fact that the definitions of the various subsets of needed to define the regions where spins 
are frozen to "— " or "-|-" ([/, J^, J^, etc) have to be adapted in the obvious way due to the lack 
of discrete-rotation symmetry of the general initial droplet T>. We will therefore skip altogether 
the proof of (6.4) and we limit ourselves to indicating the only point where some (minor) care 
has to be taken. 

The definition (4.45) of the set U is replaced by U := T>^~'^\ cf. (2.8). Let si be the vertical 
segment obtained moving downwards from the "north pole of U" until the point c where si 
meets S2, the horizontal segment obtained moving to the left from the "east pole" of U until 
c is reached. To prove the analog of (4.48), mimicking the proof given in Section 4.4.4, one 
would like to apply (6.2) in order to freeze to "— " all the spins along the two rescaled segments 
Lsi,Ls2- This is however not allowed in general, because nothing guarantees that they are 
entirely contained in LU, i.e., that c £ U (this problem does not occur for the invariant shape 
where c is the origin). The solution however is simple: one just freezes to "— " all the 
spins along the portions of Lsi,Ls2 which are inside LU, and along the shorter portion of LdU 
which connects them (call F this portion). The point is that in this situation the +/— interface 
between north and east poles follows again the corner dynamics and Theorem 3.2 is applicable. 
The freezing of "— " spins along F is equivalent to putting a hard-wall constraint in the corner 
dynamics (the interface is not allowed to cross a zig-zag path which approximates F) but this 
is irrelevant: since F is at distance of order L away from the linear profile the corner dynamics 
approaches for long times, the probability that the interface even feels the hard- wall constraint 
within the diffusive times of order we are interested in goes to zero with L (this again can 
be seen via Theorem 3.2). Other than that, the proof of (6.4) is identical to that in the V = ^ 
case. 
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7. Proof of Theorem 3.2: scaling limit for SSEP 

The first step is to discretize (3.6), so that instead of working with (f){-, •) we get $(•, •) solution 
of the analogous discrete Cauchy problem: 



(7.1) 



for every t > and x€{l,...,L — 1}. Here A is the discrete Laplacian operator: 

(A/)(x) :=/(x + l) + /(x-l)-2/(x) VxG {1,...,L-1}. (7.2) 

Note that $(x,t) = E [hxit)], with {h{t))t>o the process with generator (3.3), and that $(0,t) — 
/io(i) = ^{L,t) — hL{t) = 0. It is a standard result that solution of the discrete space heat- 
equation, converges to (p in all reasonable norms when L — )• cx) in the diffusive limit. We record 
this result here: 









= K 


$^(L,t) 




^>^(x,0) 





Lemma 7.1. 



lim max max -|$(l xLj , tL^) - L(/)(x, t)| = 0. (7.3) 

L^ootG[0,T]xG[0,l] L 



Using Lemma 7.1, we are reduced to prove 

^ 1. (7.4) 



lim ] 



max max \hx{t) — <I>(a;,t)| < eL 

te[0,TL2]s6{l,...,L-l} 



Both h.{t) and ^(^t) are 1-Lipschitz functions (for all t) so that \h.{t) — $(-,t)| is 2-Lipschitz 
and 

I max \hx{t)-^{x,t)\^ cX ^[s^ [K{t)-^{x,t)]^ ^ a^/-A . (7.5) 
ixe{i,...,L-i} J J 

As a consequence, (7.4) is equivalent to prove the following L2 convergence statement: 

Proposition 7.2. The following convergence in probability holds: 

1 ^"^ 

lim sup -^y^[hx{t)-^{x,t)f = 0. (7.6) 
■f'^°°ie[0,L2T] ^ ^ 

Proof of Proposition 7.2. The restriction of the operator A to 

= {g : {0, . . . , L} ^ M, g{0) = g{L) = 0} 

is self-adjoint (for the canonical scalar product on R-^"^ denoted in the sequel by (-, •)) and the 
family of functions 

/fc:{0,...,L}9x^y|sin(^), A: = 1,...,L-1 (7.7) 
forms an orthonormal basis of A^ of A-eigenfunctions, with respective eigenvalues 

- Afc := 2cos f ^ j - 2 < 0. (7.8) 
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As the function x i— )• hx{t) — ^{x,t) is in it can be decomposed on this basis. We use the 

(7.9) 



notation Hj^ for its A;— th coordinate (multiphed by yLj^ for convenience): 

L 



^i':=E[/^.W-^(^,t)]sin(^) 
x=0 ^ ^ 



c\2 



The quantity one wants to estimate in (7.6) is equal to 

^i;W)'<^i; sup (//f) 

We control the right-hand side by controlling each separately. 
Lemma 7.3. For every L, k £ {1, . . . , L — 1} and t > one has deterministically 

\H^\ < AL^/k. 

Moreover for any given T, w.h.p. 

\H'1\<L^^^ for every k<{\ogLf/^ and every t<L^T. 
Proof of Lemma 7. 3. The first point is easy. Using summation by parts 

T T , ^ 

kiry 



(7.10) 



y=x 



L 



Then one can check that for every x and k 

\[K{t) - $(x,t)] - [hx-i{t) - $(x - l,t)]| < 2, 



)/=a; 



sm 



kTTy 

IT 



< 



2L 

T 



(7.11) 
(7.12) 

(7.13) 

(7.14) 
(7.15) 



so that (7.11) follows. 

For the second point, first, one notices that for all k G {1, . . . , L — 1}, the functions 

L 



nkx 



L 



hh - ho 



L 



(7.16) 



Fk : ^Ml,Nl 9 /i H> ^ sin 

are eigenfunctions of C with respective eigenvalues — A^. Indeed F^ is just a linear combination 
of the coordinate function : h ^ hx (plus a constant), and it can be seen from the very 
definition of the generator C (3.3) 



using the notation 



C{Ax){h) = \{/\h){x) = \{^h){x) 
^^^^x so that (note that h e Al) 



2CFk{h) = 7I72(/fc,A/i) = ./L/2{Afk,h) 
As a consequence one can rewrite 

L 



-X,^/L/2{fk,h) 



H 



x=0 



sm 



k-KX 

~17 



hx(t) 



L 
x=0 



sm 



klTX 

~17 



hx{^) 



(7.17) 

-\kFk{h). (7.18) 
(7.19) 



and notice that := e^'^^/^H^ is a martingale. Therefore one can get the result by computing 



the second moment of Mj^ and using Doob's inequality. 
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It is not difficult to bound the quadratic variation of . Notice that 

ft 



E 



(M, 



k\2 



and that 



k=i ^ ^ 



fk7rx\ {A{h{t)){x)f 



ds < Le^^'ds 



(7.20) 



(7.21) 



so that E [(M/=)2] < L e^^Ms. Therefore (using Afc = tt'^ k"^ / L'^ {I + o(l)) uniformly for ah 
k ^ (log L) 1/3), 



sup \H^\ > a 
te[0,L2T] 



< 



sup \M^\ > a 

ie[0,L2T] 



< C 



Using this inequality for a := L'^/^ and all k ^ (logL)^/^ one gets that 



3t G [0,L2r], 3A: ^ (logL)!/^^ \H^\ > 



< 



E 



C 



fc^(logL)i/-i 



(7.22) 



(7.23) 



One can check that the right-hand side above tends to zero when L goes to infinity, which 
finishes the proof of Lemma 7.3. □ 

We now turn to (7.10): 



sup {H^f<^, sup + 



k- 



k ^ (logL)i/3 ielO.i^r] fc=|-(logL)i/3] 

The second term tends to zero (it is roug hly (logL)~i/3). The first one is less than 

k (logL)l/3 



2l_ Sr^ r7/2 / logL 



L 



This achieves the proof of Proposition 7.2 and thus also the one of Theorem 3.2. 



(7.24) 

(7.25) 

□ 



Appendix A. Proof of Theorem 3.4: scaling limit for the zero-range process 
This section follows quite closely computations in Appendix A of [24]. 

A.l. Particle system and monotonicity. For x € {—L, . . . , L} we denote rjx := hx+i — hx 

the discrete gradient of h in x. A configuration h G 0,^ can be alternatively given hj r] £ Q]^ := 
{t] : {—L, . . . , L} — )• Z}. It turns out that the zero-range process description of the dynamics 
(cf. Section 3.2) is easier to work with. 

For a more formal description of the dynamics we write explicitly its generator. For rj € Ql 

and X G {— L,... , L — 1}, we define the configuration as 



f] {x) := r]x - sg{7]x), 
?? (2; + 1) := r]x+i +sg{r]x), 

^ (.y)-=^y^ yy^{x,x + l}. 



(A.l) 
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We define V analogously for x G {— L + 1, . . . , L} replacing x + 1 in the second and third lines 
by rr — 1. The sign function sg is given by 

'l if o > 0, 



Note that the dynamics conserves the sum of the Ty's, i.e. the value of hi+i- 

Before going to the core of the proof, we need to change slightly the initial condition. In order 
to compare with the original one, one needs the following monotonicity statement: 

Proposition A.l. [Coupling] 

(i) There is a canonical way of constructing simultaneously the dynamics with generator 
(3.8) from all possible initial configurations . It satisfies the following monotonicity 
property: given and hP with h^ > Iv^ for all x, the dynamics h and h starting from h^ 
and hP respectively satisfy hx{t) > hx{t) for every t and x. Moreover, the dynamics 
started from h^ + a, a £ 7j, (a vertically translated version of h^ , including the boundary 
conditions /iq and h^^i), is simply {h{t) + a)t>o. 

(ii) There is a canonical way of constructing the dynamics with generator (A. 3) from all 
possible initial configurations rp . It satisfies the following monotonicity property: given 
rf and fp with rp. > fp. for all x, the dynamics rj and fj starting from rp and fp respectively 
satisfy r]x{t) > fjxit) for every t and x. 

Proof. The idea of the proof is using a canonical construction of the process, similarly to what 
is done in Section 2.3. It is quite classic but we perfom it here for the sake of completeness. 

• For X G {—L + l,L} we define {Tn,x)n>o and {T^x)n>o to be two IID clock processes, 
with tq^x = and t^+i^j; — Tn,x HD exponential variables of mean 2. 

• The process h{-) is constant in time except at the of the ringing times of the clock 
processes. At time Tn,x only hx is modified, as follows: hx{Tn,x + 0"'") = hx{Tn,x) + 
sg{hx-i{Tn,x) — hxijn^x))-, the other coordinates being left unchanged. At time r^j. only 
hx is modified, as follows: /ix(r^_^ + 0+) = hxir'^^x) + ^^Q^x+i{t'^^x) ~ ^^^('^n.a:))) the other 
coordinates being left unchanged. 

The reader can check that this allows to couple the dynamics from all possible initial conditions 
and that our coupling has the desired properties. This coupling induces a coupling on 77 that 
also has the right properties. 



A. 2. Changing the initial condition. We prove (3.14) working with an initial condition 
which is not the one, h^ , described in (3.11), which is random and for which the number of 
particle at a site is given by a geometric variable. The reason for this change of initial condition 
will appear in the proof of [Hi) in Lemma A. 3. We explain in this section why this implies the 
result starting from h^ . 

Given a continuous function (fP : [—1,1] — )• M with = and with a finite number of 

changes of monotonicity, set {fix)x^{-L,...,L} to be a family of independent variables with the 
following distribution: if (j)^{{x + 1)/L) — (p^{x/L) > then fjx is a geometric variable of mean 



sg(a) := < -1 if a < 0, 
if o = 0. 

V 

The generator of the chain seen in the state-space 0l is given by 



(A.2) 




(A.3) 



□ 
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L{(lf{{x + 1) / L) — 4P {x/ L)) and if {((P{{x + l)/L) — ((P{x/L)) < then —rjx is a geometric variable 
of mean L{(j)^ {x / L) — (fp{{x + 1)/L)) (with the convention that 4>^{1 + L) = 0). One sets 

x-l 
y=-L 

Note that for every e > 0, w.h.p, 

hi - <hl<hl + L^/2+£ foj. g^gj.y ^ g |_^^ . . . , L + 1}. (A.5) 

Let {h{t))t>o, {h(t))t>o be the dynamics with generator (3.8) started with initial condition /i", 
respectively, constructed using the canonical way of Proposition A.l (i). Then with high 
probability, for every t > and xE{— L,...,L} 

hx{t) - L'/^+' < K{t) < h.,{t) + (A.6) 

Therefore in order to prove (3.14) for /i(-), it is sufficient to prove it for h{-). We let fix{t) = 
hx+i{t) — hx{t) denote the gradients of h. 

A. 3. Proof of an L2 statement. For {h{t))t>Q defined above one has 
Proposition A. 2. For any t > 

-3 mx,LH)-hx{LH)f 



lim E 



(A.7) 



This result does not directly imply (3.14) (h may have a priori unbounded gradients), but it 
is not to difficult conclude from Proposition A. 2, see Section A. 4. In the rest of the section, for 
lightness of notation we write h, rj instead of h, fj. 

Before starting the proof we need some technical statements. First note, recalling the defini- 
tion of the generator (3.8), that for every x G { — L + 1, . . . , L} 

2dt^ [hx{t)] = E [sg(r/,(t)) - sg(?7.-i(i))] , 

29fE [hl{t)\ = E [2/i,(t)(sg(r/,(t)) - sg(??,_i(t)) + (|sg(r?,(t))| + |sg(r/,_i(t))|)] . 
Now some remarks: 

Lemma A. 3. The following holds (recall notations in (3.13)j.' 

(i) maxa, |(7a;(t)| is a non-increasing function oft. As a consequence 

yt > 0, Vx G {-L, . . . , L}, \qx{t)\ < \\dx4P\\oo. (A.9) 

(ii) maxj, \a{qx+i{t)) — a{qx{t))\ is a non- increasing function oft (recall cr{u) = n/(l + \u\) ). 
Then, for some C(/) = C(||9^(^0||^^ \\dl(tP\\oo) < 00 

Vt >0, VxG {-L,...,L}, |g,+i(t) - <?,(t)| <C(/)/L. (A.IO) 

(iii) For any t, the random vectors {r]x{t))x^{-L,...,L} ^"-^ (~^2:(*))x6{-l,...,l} are stochastically 
dominated by 2L + 1 IID geometric variables with mean \\dx4>^\\oo- 

Proof. For (i) it is sufficient to show that Q{t) = m.ayixqx{t) is non-increasing (by a similar 
argument one shows that Ta.m.qx{t) is non-decreasing). As the maximum over finitely many 
differentiable functions, xn&yix Ixit) possesses a right and a left-derivative everywhere and the 
right-derivative is equal to 

dtQ{t)= max dtqx{t). (A.ll) 

xSargmax q. (t) 
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For any x in max^.gargmaxq.(i), one has 

2dtq^{t) = + - 2a{q^{t)) < 0, (A.12) 

(as cy{qx{t)) is maximal), and therefore Q{t) is decreasing. 

For [ii): Using the same argument as for the point (z), we just have to note that for any fixed 
time T and xq where uia.'x.x[(T{qx+i) — a{qx)]{T) is attained one has 

2 [dt{a{qxo+i) - o-feo)}] [T) = o-'(g^„+i(T)) [a{qxo+2{T)) - o-(g^o+i(r))] 

+ a'{qxo{T)) [cr{qxo{T)) - a{qxo-i{T))] 
- [a'{qx,+i{T)) + a\qx,{T))\ [a{qx,+i{T)) - cTfe„(r))] < 0. (A.13) 

For [in): One has that L {4>^{^^^^) — 4>^{ j)) — ll'^x'^'^Hoo so that the initial configuration is 
stochastically dominated by fj^ the configuration given by 2L + 1 IID geometric variables with 
mean ||9xi;^'^||oo- According to Proposition A.l (ii), one can couple the two dynamics t] and fj 
starting from r/^ and fj^ so that r]{t) < r/(t) for all t > 0. For fixed t the law of fj{t) is the same 
as the one of fj^ as this distribution is stationary for the dynamics. The other domination is 
proved in the same way. □ 

Proof of Proposition A.2. We estimate the difference between E -p- X]xi-L(^(^' -^^*) ~ ^x(-^^*))^ 
and the same quantity at time zero, by considering it as the integral of the derivative. 



E 



-1 Y.{^[x,LH)-hx[LH)S' 



x=—L 



E 



-^3 Y.{^x,Q)-hxm' 



x=—L 



-1 / dsE[mx,s)-hx{s)f]ds 



1 ^ 

7? ^ 



x=-L+l 



x=-L+l 

E {($(x, s) - hx{s)) (a((?,(s)) - a{qx-i{s))) 

- <^{x,s){sg{r]x{s)) - sg{r]x-i{s))) 
+hx{s){sg{r]x{s)) - sg{r]x-i{s)) + ^ (|sg(r/^.(s))| + |sg(7?^._i(s)|) I ds 



13 Yl 



x=—L 



LH 



E 



qx{s)(y{qx{s)) + r]x{s)(T{qx{s)) + qx{s)sg{r]x{s)) 



Jo 



E 



{\Vx{s)\ - \sg{'r]x{s)\) 



1 



ds 



hL+i{s){sg{riL{s)) - (j{qL{s))) + -(|sg(7?_L(s)| + |sg(r?L(s)|) 



ds. (A.14) 



The second equality is obtained by expanding the product and using (A. 8) and (3.12) to estimate 
all the derivated terms. The third equality is obtained via summation by parts, it gives a term 
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that is due to boundary effect (the second one) which can be bounded as follows 



L- 



LH 



E 



/iL+i (s)(sg (r/L(s)) - (j{qL{s))) + -(|sg(r?„z,(s))| + |sg(r/L(s))|) 



ds 



< CL^\l + E|/iL+i|) = 0{L-^'^). (A.15) 

Indeed hL^i{t) is constant through time and is the sum of 2L + 1 independent variables. The 
mean of this sum is and the variance of each term is bounded as we supposed cfP to be smooth. 
The variance of is thus 0{L). We can also neglect the second term in the first line as 



E 



L+l 1 L+1 

x=—L J x=—L 



(A.16) 



where the last equality is easy to obtain once noticed that is the sum of (L + x) independent 
geometric variables with bounded variance. 



Set 



A{x,s) 



-qx{s)a-{qx{s)) + r]a:{s)a{qa:{s)) + gx(s)sg(r/x(s)) - {\Vx{s)\ - |sg(77^(s)|). (A. 17) 



From the previous equations one gets that 



E 



1 ^ 

7? E 



i'^{x,Lh)-h,{Lh)f 



-L+l 



1 

L3 



L 

E 

=-L+l 



E[A{s,x)]ds + o{l) 



(A.18) 



To understand better the rest of the proof, the reader should notice that if i'r]x{s))xe[-L,...,L} 
were distributed like geometric variables it would be possible to factorize E [A{x, s)] in a product 
of negative sign and from equation (A.18) the proof would be over. Indeed, for q > and rj 
distributed like a geometric variable of mean u > (or —i] is distributed like a geometric variable 
of mean —u > 0), 



E [-9^(9) + v(^{q) + qsgiii) - {\v\ - \sg{v)\)] = -{q- u)Hq) - (^(u)) < o. 



(A.19) 



(recall that ct(-) is an increasing function). It is not true in general that r]x{s) are geometrically 
distributed for s > but this is reasonable to think that their distribution is close to geometric: 
as the system mixes locally in finite time, what one should observe on finite but large windows 
is close to an equilibrium measure, and from [1] it is known that the only (infinite-volume 
translation invariant) equilibrium measures for the zero-range process are convex combinations 
of products of geometric variables. Most of our efforts will therefore be focused on proving 
convergence to the infinite volume measure for a space-time averaged version of the probability 
distribution of the r]x{s) (using this space-time average is somehow crucial for the proof to work). 

As the limiting object is an infinite volume measure, it is somehow more convenient to consider 
r]{s) as an element of by periodizing it: for the system of size {2L-\-l) one sets 

^x+k{2L+l) ^x 

for every € Z,x G {—L, . . . ,L}. For y G Z one defines 6y to be the shift operator r] i— )■ 6xT] 
defined by 

Vx G Z, {6yri)x := Vx+y (A.20) 

We define for each L > the measure /U^ on Z^ our space-time averaged measures by its 
action on local functions (for G N we call f{r]) a A'-local function if / is bounded and can be 
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written as a function of 'r]^[-K,K]i / is a local function if there exists a K such that it is i^-local): 



1 1 
tl2 2L + 1 Jo 



J2 f{ey{v{s)))ds 



y=-L 



(A.21) 



We want to prove that any limit point (when L — t- oo) of /xf is an equilibrium and use this 
information to bound the right-hand side of (A. 18). 

We introduce some notation to describe the limiting measure. For n G M define to be a 
measure on 77 = {r]x)xez such that the rjx are IID geometric variables of mean n if u > while 
the —r]x are IID geometric variables of mean —u ii u < 0. If is a probability measure on M 
define 

p'^ ■= [ p^i,{du). (A.22) 



Proposition A. 4. Fix t > 0. For any subsequence of {fi^")n>0! it is possible to find a sub- 
subsequence (//^ ")n>o that converges locally to with v a probability measure on M with support 
included in [—\\dx(i)^\\ooi ||<9x<^'^||oo]; in the sense that for any local function f 

hm pf"(/) = p-(/). (A.23) 

As a consequence for any local function f 

limsup^,^(/)< max /,"(/). (A.24) 

L-¥oo u£l-\\dx<p'-'\\rx,,\\Ox4> Wool 

Remark A. 5. Note that the convergence does not hold in the topology induced by the total 
variation distance: indeed fij^ give mass one to L-periodic rj whereas these configurations have 
mass zero for the limiting measure. 

Proof of Proposition A. 4- For any fixed K > 0, the sequence of laws of {r]x)xe[-K,K] under pf" 
is tight by Lemma A. 3 (Hi) and hence we can extract a converging subsequence. By diagonal 
extraction it is possible to extract a subsequence of L.„ and a family of measures (^/<)/< ^ on 
such that the law of {'nx)xe[-K,K] under p^ " converges to fi^ for all K. By construction 
for H ^ K, ph projected on is equal to fix and by Kolmogorov extension theorem, 

there exists a measure fi on such that fi projected on Z^"^'^! equals hk for all K. One has 
therefore for all local function / 

lim pf"{f) = p{f). (A.25) 

n— >oo 

We have to show that p can be written as p'^. First one remarks that p^" is translation invariant, 
so that p is too. A second point to make is that p -almost surely all the r]x (that are not equal 
to zero) have the same sign. Indeed 

p (3x, X G Z, •^^x'^lx < 0) = lim p (3x, x! G [— -f^, -f^], ^x^x < O) 



lim lim /if" (3x,x' G [-K, K],r]xr]x < O) (A. 26) 



and 

p\ (3x,x'g [-K,i^],77.r/; <0) 
1 



tL^{2L + l) 



E 

y=-L 



[3x,x' G {-K + y,K + ylT]x{s)r]x'[s) < O] ds. (A.27) 
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One realizes easily that 

L 

'^{3x,x'(i{-K+y,K+y],n:,'r]^,<iA\ (A.28) 

y=^L 

is upper bounded by (2i^ + 1) times the number of changes of sign in {r]r^x^<y_i^i^j^^. Prom 
the definition of the dynamics, a transition can only lower the number of changes of sign. Its 
initial value is smaller than the number of changes of monotonicity of ^ (which is assumed to 
be finite) plus one (the "plus one" can come from periodizing). Therefore 

L 

^ P[3x,x' G [-ir + y,i^ + y],77,(s)r?,,(s) <0] < 2JfC(/). (A.29) 

y=-L 

A third point is to show is that [i is an invariant measure for the infinite volume dynamics 
(the infinite volume version of (A. 3), call its generator £°°). For / a i^- local function one has 
(for L> K large enough) 

^i{C°^f) = ^^j^ f (^"(^ ° Oy)ivm (A.30) 

y=-L 

For y G [—L + K,L — K] the infinite volume generator applied to / has the same effect as the 
finite volume generator so that 

f E [C^if o ey)iv{s))] ds = f dsE [(/ o 6y)Hs))] ds = E [(/ o ey){7^{t)) - (/ o 0,)(r?(O))] . 

JQ Jo 

(A.31) 

Therefore 



y=-L+K 
.LH f-L+K-1 



^W^Wnf £ + E IE(/:°°(/o0,)(r?(s)))d. = O(l/L). (A.32) 

^ JO \ y=-L y=L-K + lJ 

As a consequence, for any local function 

m(/:°°/)= hm ^f"(£°°/) = 0. (A.33) 

n— >oo 

Restricted on the event r]x have all the same sign, is the generator of the zero-range process 
with one type of particle and therefore is a translation invariant measure for the zero-range 
process. From [1, Theorem 1.9] one can write ji = p'^ for some v. By Lemma A. 3 {in) — (iv), 
under at time zero r] is dominated by a IID family of geometric variables of mean ||9a;0'^||oo 
and so is —r]. This implies the claim on the support of i/. 

The second point of Proposition A. 4 is standard; we include its proof for completeness. Given 
a local / one can extract a subsequence Ln such that 

lim /i^(/) = limsup;u,^(/). (A.34) 

From Ln one can extract a subsequence L'^ such that /x^ " converges to p"^ so that 

lim /.^(/) = lim /if"(/) = / p'^ifHdu), (A.35) 
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which ends the proof. 

Fix / to be a large fixed integer. Set for y G Z 

By ■.= {l + y,...,l + y}. 



□ 



(A.36) 



For notational convenience, similarly to t] with (A. 20), one considers now periodized version 
{qx{s))sez of q{s) and {A{x,s))x(zz of A{-,s). 

Now, one uses Proposition A. 4 to control each term in ¥,^A{s,x). 



Lemma A. 6. 



g^limsup -^E / XI ( 7 XI ^^■(•5)sgfe(s)) j - qyis)(T ij^Vxis)] ds 
° y=-L \ x£By J \ xeBy / 

Proof of Lemma A. 6. Fix I > 0. For L large enough, any all y G {—L, . . . ,L — I}, 

J qx{s)sg{i]x{s)) j - qy{s)a Ij ^ ^^(^) 1 
, ^6Ba / \ x&By J 

( 7 X sgfe(s)) I - ^ I 7 X 1 

V XeBy J \ X&By J 



-- 0. 

(A.37) 



+ msix\qx{s) - qy{s)\. (A. 38) 

X&By 



Moreover, uniformly in y € {—L, . . . , L — /}, as a consequence of Lemma A. 3 (ii) 

\qxis)-qyis)\= 0(1/1). (A.39) 



max 

y&{-L,...,L-l},xeBy,s>0 ' 



The contribution of y G {L — Z + 1, L} to the sum under the integral in (A.37) is 0{l). Therefore 
summing over y G { — L, . . . , L}, integrating over s and taking expectation one gets 



E 



X] I 7 X] (lxis)sg{Vx{s)) I - qy{s)a ( y ^ '?^(' 



y=-L \ X&By 



xeBy 



< {max\qy{s)\) 



y 



E 







X 7 X sgfe(s)) - ^ 7 X ^-(' 



y=-L \ x&By 



X&By 



{max\qy{s)\)tL^{2L + 



+ 0{lL^ 



+ 0(/l2) (A.40) 



where //f is defined in (A. 21). Therefore, the proof of our statement is finished provided one 
proves 



lim lim sup //^ 



\ x£Bo I \ xeBa j 



0. 



(A.41) 
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From Proposition A. 4 one has 



lim sup /ij 

L— >oo 



x&Bo J \ xeBo 



< sup 

o<-u<||a^(AO 




(A.42) 



and one can check that the right-hand side term tends to zero when I tends to infinity: we note 
that under p", for every x one has p^{sg{i]x)) = cr{u), and the law of large numbers tell us that 

the two terms j YlxeBy ^s(^x) and a (^j ^x<^By Vx^ have the same limit when / tends to infinity. 
However, because of the sup over u one needs more quantitative estimates than the law of large 
numbers to conclude. For instance we can get them by the use of second moment method; we 
leave the details to the reader. □ 



Similarly to Lemma A. 6 one shows that 
Lemma A. 7. 



1 r^^* ^ f(2T -\- 1 

lim limsup-^ / V E(G(7?(s))) = lim lim sup ^ ' f4{G{r])) = (A.43) 



y=-L 



where 



Gir,) 



1 



x&By 




(A.44) 



Proof. The proof is very similar to that of Proposition A. 6, the only additional technical point 
being that the function G{r]) is not bounded so that one cannot use directly Proposition A. 4. 
However stochastic domination given by Lemma A. 3 {in) allows us to get the same conclu- 
sion by considering the function r] i— )■ mm{G{r]), K), and letting K tend to infinity afterwards. 
Altogether one gets 



limsup^t'(G) < sup p"(G). 



L— >oo 



0<u<\\d^4>o\\ 



(A.45) 



We end the proof in the same way that for the previous Lemma remarking that 



p'^ilVxl - \sgivx)\) = ua{u). 



□ 
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Now we are ready to conclude: 

L 



y=-L 
L 



\ -Qyis)cr{qyis)) + j i Vxis)a{q:r{s)) + qx{s)sg{ri:^{s)) - (|r?^(s)| - |sg(r/^.(s)|) 
\xeBy 

y=-L \ x&By J \ x&By ) 



R{s,l,L)- Y 

V=~L 



\ Y ^-(^) U 7 

xeBy I \ 
Qy{s) - ij Y 1 



Y Vx{s)j +R{s,l,L) 

(^ilyis)) - o- J Y 'nx{s) 



XGBy 



(A.46) 



where 



R{s,l,L) = - Y 

y=-L \x&By J 

I 7 Y 9^(s)sgfe(s)) I - Qy{s)o- ( J Y '?^(' 

\ ^^By J \ X&By 

I 7 X] ~ |sgfe(s)| I - I 7 X] 1 ^ I 7 2Z 

\ X£By J \ X&By J \ XdBy ^ 



(A.47) 



and the second term is non positive {a — b and cr{a) — a{b) have the same sign). 

According to (ii) — {iii) in Lemma A. 3 (to control the first term) and Lemmata A. 6 and A. 7 



1 z-^'* 

limlimsup— ^ / K R{s,l, L)ds = 0. 
L-i>0 Jo 



This implies 



1 f 

limsup— ^ / > 'EA(x,s)ds<0 

L^oc L'^ Jo _r 



y=-L 



and therefore the result that we want to prove, from (A. 18). 



(A.48) 



(A.49) 



□ 



A. 4. Concluding the proof of Theorem 3.4. It is not hard to transform the L2 statement 
of Proposition A. 2 into the desired "almost sure" statement: 



Proposition A. 8. For any e > 0, t > 0, w.h.p 

1 



max ^\^{x,L'^t) - hx{L^t)\ < e. 
ve{-L,...,L+i} L 



(A.50) 
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Proof. Also here we write h for h. Note that from Lemma A. 3 {in) the random vector {\r]x{t)\)x^{-L,. 
is stochastically dominated for every t by a vector of IID time-independent geometric variables. 
This implies that there exists a constant C such that for any t > 0, w.h.p. 

\hx(t) — hy{t)\ < C\x — y\ for every X, y G {—L, L}, jx — y| > log L (A. 51) 

(this can be proved by using large deviation estimates and a union bound onx,y £ {—L, . . . , L}). 
Moreover Lemma A. 3 (i) ensures that is always Lipschitz so that (A. 51) holds also for 

<^{;t)-h.{t). 

With (A. 51) and L large enough, one has 

max \^(x, L^t) - hJL^t)\ > eL) (A.52) 

a;e{-L,...,L+l} J 

mx,Lh)-hx{Lh)\^>'-^\ (A.53) 

^a;e{-L,...,L+l} J 

SO that the left-hand side event has small probability when L is large, otherwise Proposition A. 2 
would be false. □ 

A. 5. Laplacian bounds. Recall that <I>(a;,t) is the solution of the Cauchy problem (3.12). We 
want to bound $(x,t) above and below with the solution of a suitable heat equation. For this, 
we will suppose that the function (fp , through which the initial condition $o for $(x, t) is defined, 
is concave on [—1, 1] (in addition to the assumptions required for Theorem 3.4). One defines the 
evolution $i(x,t) as the solution of 

'dt^i{x,t) =\/^^i{x,t) 

$i(-L,t) =$(L + l,t) = (A.54) 
$i(x,0) =$o(x) 

for t > 0, x S {— L + 1, L}. Also we define ^2{.Xi t) as the solution of the analogous equation (with 
the same boundary values) where the discrete Laplacian is multiplied by {l/2)a' {\\dx(j)^\ 
1/(1 + \\d. 



-'xH^ Woo) 



Proposition A. 9. For every t > every x £ {—L, . . . , L + 1} one has 

^i{x,t) <<^>{x,t) < ^2{x,t). (A. 55) 

Proof. We prove the upper bound, the lower one being very similar. Suppose that the result 
does not hold and set 

T := max{t | ^{x, t) < $2(2;, t) for every t < T, x G {-L, ...,L + 1}}. (A. 56) 

Note that by property of the heat-equation, <I>2(x,i) is a strictly concave function of x for all 
positive t (except in the case where one starts from the flat initial condition but in that case the 
statement is trivial). Let xq be such that 

^{xo,T) = ^2{xo,T). (A.57) 

Then one remarks that q^giT) — qxo-i{T) < (by strict concavity of <I>2(";^)) and that by 
Lemma A. 3 max^. |(?x(i)l — l|f^x''/''^l|oo so that 

a{qxo{T)) - a{qxo^i{T)) < (g,„(r) - qx,^i{T))a' {\\dxM\oo) (A.58) 
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(since cr'(-) is decreasing on R"*") and hence 

2dt[^2 - m^o,T) = a'{\\d,M\oo)^Mx,t) - a{q,,{T)) + a{q^^^i{T)) 

> <j'i\\d.M\oo)[iMx + l,t) + $2(x - 1, t)) - + 1, t) + - 1, t))]. (A.59) 

Since the last expression is non-negative, one has < ^2ix,t) on an interval [T, T + £{x)] 

for some £{x) > 0, for every x G {—L, . . . ,L + 1} and that concludes the proof since the only 
possibility is that T = oo. □ 
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